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1.  INTRODUCTION 


In  this  report,  we  study  aspects  of  the  role  of  plasma  turbulence  in 
late-time  striation  behavior  and  satellite  communication  effects. 

Sections  2  through  4  are  general  and  conceptual.  We  begin  in  Section  2 
by  deriving  a  basic  picture  of  particle  diffusion  due  to  quasi-neutral  plasma 
modes.  In  Section  3,  we  relate  the  form  of  the  diffusion  coefficient  to 
linear  and  non-linear  plasma  theory.  The  results  for  the  form  of  the  diffu¬ 
sion  coefficient  are  shown  in  Eqs.  (4),  (7)  and  (8).  In  Section  4,  we  discuss 
as  it  might  apply  to  E  x  B  modes  in  the  absence  of  other  modes  gaining  insight 
into  some  of  the  limitations  of  the  concept,  as  well  as  possibly  into  bifurca¬ 
tion. 

In  Sections  5  through  8,  aspects  specifically  related  to  drift-dissipa¬ 
tive  modes  are  presented.  In  Section  5,  we  discuss  diffusion  due  to  drift- 
dissipative  modes  (DDM)  showing  that  the  contribution  at  a  given  altitude 
along  a  field  line  should  not  be  shorted  out  by  the  behavior  at  other  altitudes, 
and  that  DDM  diffusion  should  act  to  limit  striation  bifurcation.  Perpendicular 
diffusion  coefficients  and  rates  of  perpendicular  diffusion  due  to  DDM  along 
in  the  absence  of  background  convection  and  other  diffusive  mechanisms  are 
presented  as  functions  of  background  parameters.  It  is  shown  here  (and  in 
Appendix  B)  that  effects  due  to  an  electric  field,  E^,  in  the  direction  of  the 
background  striation  density  gradient  and  perpendicular  to  B  limit  the  growth 
and  diffusion  as  long  as  the  velocity  cEx/B  is  larger  than  an  appropriately 
averaged  ion  diamagnetic  drift  velocity  =  (v*^)ay.  In  Section  6,  diffusion 
coefficients  from  drift-dissipative  modes  are  characterized  for  use  in  computer 
simulations.  In  Section  7,  observational  aspects  of  drift  waves  in  the  ionos¬ 
phere  are  considered,  and  in  Section  8,  specific  experimental  situations  are 
treated.  These  latter  include  aspects  of  the  BUARO  release,  STRESS,  and 
three-meter  backscatter. 

In  Section  9,  conclusions  and  experimental  suggestions  are  presented. 
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2.  BASIC  PICTURE  OF  TURBULENT  DIFFUSION 
FOR  QUASI -NEUTRAL  MODES 


The  basic  picture  of  turbulent  diffusion  for  quasi -neutral  modes  in  a 
plasma  can  be  gotten  by  assuming  that  both  ions  and  electrons  behave  according 
to  the  Vlasov  equation  with  particle  collisions  neglected.  This  would  appear 
to  be  valid  as  Tong  as  the  diffusion  due  to  turbulence  dominates  the  diffusion 
due  to  ordinary  collisions. 

We  consider  modes  which  vary  on  a  relatively  fast  time  scale  compared  to 
the  turbulent  diffusion.  The  magnetic  field  B  is  chosen  to  be  in  the 
z-direction  and  the  background  density  gradient  is  in  the  x-direction.  Then 
the  Vlasov  equation  for  each  species  can  be  written 


3f  + 
3T  +  V 


if  A  qs  ♦  v  *  3f  .  qs  ±  3f  _  _ 

*ar  +  m  v  ®*3v  +  m  3v  "  0 
s  s 


with  f  the  particle  distribution  function  in  r,  v  and  t  space;  q$  =  species 
charge,  mg  =  species  mass;  E  =  electric  field;  and  subscript  e  =  electrons, 
subscript  i  ==  ions. 

On  writing 

f  =  fo  +  fl 


(with  f  the  slowly  varying  background  distribution,  and  fj  the  distribution 
associated  with  the  turbulent  modes)  and  averaging  Eq.  (1)  for  a  time  long 
compared  to  the  time  scale  for  variation  of  the  turbulent  modes  one  obtains 

3f  3f  q  +  *  3f  .  q  3f  q„  /  3f,\ 

o  .  -*•  o  ,  Ms  v  x  8 _ o  .  Ms  f  o  a.  s  /f  .  1\  _  n  r?\ 

W  +  v*3 r  +  ms  c  3v  +  mg  ^*3v  +  m$  vl  3v  /  "  0  '  ^ 

where  <  >  denotes  the  time  average  as  well  as  a  summation  over  modes.  Multi¬ 

plying  through  by  v  and  integrating  over  velocity  space  yields 


xr(nv)  +  V  (nvv  +  c 


> +/ \  * -c-  *  ^d-  +/^vE*^rd-  +/^(E!‘^r)  d.  =  0 


where  is  the  particle  random  energy  stress  tensor, 
cation  this  becomes: 


After  some  simp! if i- 


10 


■  sses 


3^(nv)  +  V-(nvv  +  <(>)-  ~r  (nv^ 


where  v  o  ^  »  vp>  yfith  p  being  the  pressure,  subscript  "1"  denoting  a  first- 
order  quantity,  n  is  the  ionized  number  density  and  v  is  now  the  fluid  veloci¬ 
ty.  If  the  background  motion  has  a  time  scale  t,  velocity  v  and  length  scale 
L,  then  one  expects 

lqs!B  »  i 

cs  me  t  » 
s 

This  allows  one  to  neglect  the  terms  in  a/at ( nv)  and  v*(nw).  Further  we 
expect  no  background  pressure  gradient  in  the  y-direction.  Then  on  taking  the 
y-component  of  Eq.  (3)  one  has: 


nvx  3  !  <Ely"l>  +  I  Eyn 


Convection  in  the  x-direction  is  associated  with  the  term  cEyn/B.  It  is  clear 
that  diffusion  in  the  direction  of  the  density  gradient  (x-direction)  must  be 
associated  with  the  term  (c/B)<E^y  n^>.  This  term  is  identical  for  ions  and 
electrons  provided  the  modes  are  quasi -neutral ;  hence  diffusion  occurs  locally 
at  a  given  point  in  three-dimensional  space  independent  of  ambi-polar  electric 
fields  or  the  short-circuiting  effects  which  can  arise  with  ordinary  particle 
collisions^  in  the  inhomogeneous  ionosphere. 
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3.  DETERMINATION  OF  THE  FORM  OF  DIFFUSION 


If  plasma  turbulence  is  in  a  quasi-stationary  state,  i.e.,  the  mode 
amplitudes  are  changing  slowly  compared  to  the  mode  frequencies,  the  product 
<E^nj>  is  determined  by  a  balance  of  linear  and  nonlinear  growth  and  damping 
processes.  In  place  of  solving  strongly  nonlinearly  coupled  equations  for 
mode  amplitudes,"  estimates  for  can  be  obtained  by  using  the  relation 

between  Ej  and  nj  for  linear  growth  (i.e.,  at  amplitudes  below  saturation) 
and  using  nj  as  estimated  at  saturation.  Computer  simulations  of  universal 
modes^’^3^  suggest  this  to  be  a  reasonable  procedure. 

It  is  conceptually  simpler  to  use  electrons  to  establish  the  linear 
relation  between  E^  and  n^  because  finite  Larmor  radius  effects  are  negligi¬ 
ble  provided  kpg  <<  1  where  k  is  the  model  wave  number  perpendicular  to  B  and 
Pg  is  the  electron  Larmor  radius.  If  this  is  the  case,  the  mean  electron 
drift  in  the  x-direction  is  cE^/B.  If  the  mode  has  a  growth  rate  Y  and  flow 
for  electrons  perpendicular  to  the  magnetic  field  is  incompressible,  then  the 
fractional  growth  per  unit  time  is 


-1  ly  an  _ 
n^  B  3x” 


Y 


(The  minus  sign  is  used  since  an  unstable  mode  will  lead  to  a  flux  of  plasma 
opposite  to  the  density  gradient.)  This  results  in 

2 

"V'WW)  +  f  V  •  (4) 

Provided  that  n^  is  randomly  varying,  this  form  for  the  diffusion  equation 

suggests  that  a  close  relationship  between  scintillation  and  diffusion  exists 

since  scintillation  in  the  weak  scattering  limit  is  a  linear  function  of  the 

(41 

three  dimensional  power  spectrum.'  '  Of  course,  scintillation  may  also  arise 

if  n,  is  not  randomly  varying. 

1  (5) 

For  strong  turbulence,'  '  with  y  comparable  to  the  wave  oscillatory 
frequency  u>,  n^  is  limited  by  the  condition  that  its  perturbed  density  gradi¬ 
ent  in  the  direction  of  the  mean  density  gradient  (the  x-direction)  becomes  of 
the  order  of  the  mean  density  gradient.  Then  one  has 
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(5) 


For  weak  turbulence,  with  y/u> «  1,  if  in  addition  there  exists  a  sur¬ 
face  of  wave  vectors  k'  in  k-space  such  that 

%  =  +  * 

(with  a),  denoting  frequency  as  a  function  of  wave  vector)  it  has  been 
estimated'  '  that  the  diffusion  coefficient  is 

Dx  *  <y2/u>k2>  . 
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4.  DIFFUSION,  E  x  B  MODES,  AND  BIFURCATION 


Application  of  the  notion  of  turbulent  diffusion  to  E  *  B  modes ^  in 
the  absence  of  other  modes  leads  to  an  appreciation  of  the  approximate  nature 
of  the  concept  especially  as  applied  to  ionospheric  scintillation  as  well  as 
to  insight  into  the  process  of  bifurcation. 

Turbulent  'diffusion  is  a  concept  based  on  the  separation  of  space  scales 
between  the  background  density  gradient  and  the  turbulent  modes  which  serve  to 
diffuse  it.  For  E  *  B  modes,  wavelengths  approaching  <(l/n)(dn/dx)>"*  are 
unstable  so  that  the  separation  of  space  scales  does  not  generally  apply. 
However  one  interesting  and  illuminating  paradox  involving  E  x  B  modes  and 
turbulent  diffusion  can  still  occur. 

We  note  that  with  E  *  B  modes  alone,  the  quantity 

f  n2  d V  , 

yall  space 

with  V  indicating  volume,  must  be  conserved. ^  To  see  this  we  observe  for 
electrons  in  the  absence  of  finite  Larmor  radius  effects. 


3n  cE  x  B 
at  b2 


Vn  =  0  , 


or  multiplying  by  n,  using  V-(E  *  B)  =0  and  integrating  over  all  space 


a 

at 


However  for  relatively  short  wavelength  E  x  B  modes  it  is  still  possible  to 
derive  a  diffusion  coefficient  since  whenever 


and 


one  has 


v  I  dn 
y  n  '3x  ’ 
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Then: 


D  ‘  "short  £  TiTij"  ^  ° 

wavelengths"  ^ 


cE  *  B 


?n  *  V‘(Dj_7n) 


Jn2  dV  =  -J D±v 


Vn  dV 


The  resolution  of  the  paradox  lies  in  noting  that  the  diffusion  equation 
only  applies  to  length  scales  greater  than  the  (short)  turbulent  mode  length 
scale.  One  expects  that  turbulence  involving  modes  at  shorter  scale  lengths 
should  act  to  decrease 


when  n  is  restricted  to  Fourier  components  at  wavelengths  greater  than  the 
turbulent  wavelengths,  precisely  because  Jn2  dV  is  rigorously  conserved  over 
all  wavelengths  for  E  *  B  turbulence. 

For  insight  into  bifurcation  we  consider  a  striation  of  thickness  L' 
with  gradient  scale  length  lx  (both  dimensions  measured  perpendicular  to  B.) 
[Typically  in  bifurcation  one  would  expect  a  striation  to  break  up  into  wave¬ 
lengths  of  2L'/3  (see  Figure  1).  Certainly  the  picture  of  turbulence  super¬ 
posed  on  a  laminar  background  of  much  larger  scale  length  is  inappropriate.] 
For  a  plasma  in  the  limit  v^n  >>  x"*,  where  v.^  is  the  ion-neutral  collision 
frequency  and  x  is  the  time  scale  for  plasma  velocity  variation,  simula¬ 
tions^  indicate  that  bifurcation  is  suppressed  for 

L  v 

-rr-^  <  400  ,  (9a) 
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with  vQ  the  neutral  wind  velocity.  For  the  equations  used  in  the  NRL  simula¬ 
tions  modal  growth  rates,  y,  can  be  estimated  by 


with  o  expected  to  be  less  than  unity.  The  result  with  a  =  1  is  based  on  an 
“ideal"  geometry  with  uniform  ion  neutral  velocity  vQ  in  the  x-direction  anti¬ 
parallel  to  the  background  density  gradient  and  perpendicular  to  the  direction 
of  spatially  oscillatory  variation  of  the  E  x  B  mode.  For  non-ideal  cases  one 
expects  o  <  1.  Further  the  eigenmode  is  assumed  to  have  no  variation  in  the 
x-direction  (local  approximation). 

For  stability  with  k  =  3*/L'  this  requires 


The  similarity  of  this  estimate  to  the  computer  results  provides  two  results. 

1.  The  Reynolds  number  criterion  L  vq/Dj_  <  400  for  bifurcation 
suppression  may  well  be  a  direct  reflection  of  the  suppression  of 
striation  instability  at  sufficiently  small  values  of  L ,  vA/0  . 

(9)  0 

Although  L^/L'  <  1  for  the  simulations'  '  one  would  expect  that  a 
more  realistic  value  of  a  (a  <  1)  would  work  to  increase  the 
critical  value  of  LjVQ/D  , 

2.  The  presence  of  the  factor  (L^/l')  on  the  right  hand  side  of  Eq. 
(9a)  suggests  that  bifurcation  in  the  numerical  simulations  could  be 
a  function  of  the  extent  of  the  grid  (in  the  y-direction) 
perpendicular  to  the  initial  density  gradient  compared  to  the 
initial  density  scale  length,  e.g.,  the  larger  the  extent  of  the 
grid  in  the  y-direction  the  larger  the  possible  values  of  l'  and  the 
smaller  the  value  of  LjV  /D^  at  which  bifurcation  is  suppressed. 
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5.  DIFFUSION  WITH  DRIFT-DISSIPATIVE  MODES 


We  have  previously  treated  some  aspects  of  drift  dissipative  modes. 

A  new  discussion  of  the  basic  physical  mechanism  involved  in  driving  the  mode 
is  presented  in  Appendix  A. 

Drift  dissipative  modes  are  quasi-neutral.  Hence  electrons  and  ions 
diffuse  identically  across  the  magnetic  field  without  an  ambipolar  electric 
field.  This  means  that  diffusion  can  take  place  at  one  altitude  along  a  mag¬ 
netic  field  line  without  "shorting  out"  or  limiting  from  contributions  at 
other  altitudes. 

The  largest  growth  rates  occur  for  k^.  *  l/10^  provided  that  ion- ion 
collisions,  v.^,  are  sufficiently  weak  that  0.3  <,  vin,  the  ion-neutral 
collision  frequency.  Further,  in  this  case,  one  expects  modes  to  be  damped 
whenever  the  wave  frequency  is  comparable  to  vin,  i.e.,  with  wave  numbers  such 
that 


<  L .  v.  /C 
X  in  s 


9 


(with  C$  the  ion-sound  speed  =  /Tg/m. ,  and  temperatures  including  Boltzmann's 
constant).  Hence  for  a  striation  of  thickness  l'  with  ky  =  3w/L'  for  bifur¬ 
cation,  short  wavelength  turbulence  is  expected  for 


L'  >>  3n 


P,C 


is 

L_LVin 


This  implies,  as  shown  in  Figure  2,  that  drift  dissipative  modes  can  act  to 
diffuse  striatlons  on  a  large  scale  or  averaged  sense. 

Since  bifurcation  is  a  large  scale  length  phenomenon  the  diffusion  due 
to  drift-dissipative  turbulence  should  be  relevant  in  limiting  it  according  to 
Eqs.  (9)  or  (10).  On  letting  vts  ■  ^%/m  ,  using  Eqs.  (7)  and  (9a)  with  y  * 
0.1  for  "vigorous"  growth  and  kx“  *  p^  the  criterion  for  bifurca¬ 
tion  becomes 
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a.  Original  density  contours  n^  and  n2 


b.  Density  contours  with  turbulent  diffusion 


c.  Averaged  density  contours,  diffusion  on  a  large  scale. 


Figure  2.  Turbulent  diffusion:  (a)  original  density  contours 
n,  and  n2;  (b)  density  contours  with  turbulent  dif¬ 
fusion;  and  (c)  averaged  density  contours,  diffusion 
on  a  large  scale. 
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P-j  V.  - 

iL>  40 .  (11) 

o 

2 

Since  characteristically  vtl./vQ  ~  10  and  Pj  s  5  x  10  cm,  bifurcation  is  lim¬ 
ited  to 

>  2  km. 

provided  "vigorous"  drift  dissipative  turbulence  is  present.  (Here  vQ  is  the 
neutral  wind  velocity.)  If  vigorous  drift  dissipative  turbulence  does  not 
break  out  until  some  scale  length  L  *  <  2.0  km  it  is  clear  that  bifurcation 
will  continue  ddwn  to  L^*. 

We  have  previously  indicated  rough  estimates  of  growth  rates  and  diffu¬ 
sion  rates  for  drift-dissipative  modes  in  the  absence  of  electric  fields  in 
the  direction  of  the  background  density  gradient.^  Modal  growth  rates  taking 
into  account  finite  ion  larmor  radius  effects  and  electron-ion  collisional 
effects  are  derived  in  Appendix  C.  The  results  are  approximately  valid  for 


ve  >  u>  . 

(Here  Lz  =  length  scale  parallel  to  3.) 

Figures  (3)  to  (7)  and  (13)  to  (16)  show  the  diffusion  coefficient  as 
given  by  the  strong  turbulence  result,  Eq.  (7),  which  is  derived  on  the  basis 
of  Eq.  (4)  and  expression  (5).  We  note  that  the  rough  estimate  for  nj 
indicated  in  expression  (5)  is  also  derivable  by  consideration  of  a  wave 

scattering  of  particles^1^  which  can  lead  to  particle  diffusion  perpendicular 

- *  2 

to  the  magnetic  field'  '  and  hence  wave  damping.  The  value  of  k^  is  taken 

to  be  (p.L  )-l  in  keeping  with  previous  results. 

1  2 
The  reciprocal  of  the  time  scale  due  to  diffusion  is  estimated  as  D^/L^. 

These  values  are  shown  in  Figures  (8)  to  (12)  and  (17)  to  (20)  as  functions  of 

background  parameters.  We  caution  that  a  given  background  density  gradient 

2 

with  diffusion  rate  D^/l^  need  not  actually  be  diffusing.  Due  to  the 
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interplay  of  convective  and  diffusive  processes  as  indicated  by  Eq.  (4)  it  is 
possible  for  a  density  gradient  with  diffusion  present  to  be  increasing  with 
time  rather  than  decreasing.  As  shown  by  barium  cloud  releases  and 
simulations,  convective  processes  alone  frequently  tend  to  increase  density 
gradients  and  hence  can  be  effectively  anti-diffusive.  Estimates  for 
gradients  at  striation  tips  have  been  derived  from  the  balancing  of  convective 
steepening  and  diffusion. 

Modal  parameters  associated  with  the  diffusion  are  shown  in  Tables  1 
6  —3 

through  5  for  n*=  10  cm  and  neutral  density  ranging  from  neutral  density  = 
10^  cm“^  (Table  1)  to  neutral  density  =  10^  cm“^  (Table  5)  in  decades. 

Values  of  Dj_  and  D^/L^  for  barium  clouds  are  shown  respectively  in 

Figures  (21)  and  (22)  and  Figures  (23)  and  (24).  For  these  figures  v_.„  = 

-10  in 
10  n  and  L  *  25  km.  There  is  no  instability  and  hence  no  anomalous 

n  1  io  _3  6-3  g  _3 

diffusion  for  n„  >  10  cm  at  n  =  10  cm  or  for  n„  >  3.16  x  10  cm  at  n 

7  -3  n  n 

»  10 '  cm  . 

A  major  limit  to  the  applicability  of  drift-dissipative  mode  turbulence 
has  been  the  uncertainty  of  effects  due  to  an  electric  field,  Ex,  in  the  di¬ 
rection  of  the  density  gradient  perpendicular  to  the  magnetic  field.  This 
should  be  a  significant  factor  in  the  occurrence  of  drift  dissipative  modes 
because  striations  if  driven  by  an  ambient  electric  field  tend  to  elongate  in 
the  direction  perpendicular  to  the  field.  Hence  along  most  of  their  perimeter 
the  electric  field  is  in  the  direction  of  the  density  gradient. 

The  problem  is  treated  in  Appendix  B.  In  the  limit  |vg|  »  (v1-*)av  with 
ve  3  cEx/B  and  (v^*)  3  -Tj/(m^uc^$0  with  a  typical  scale  length  for  the 

region  of  density  gradient,  the  growth  rate,  y,  has  an  additional  contribution 
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Diffusion  Coefficient 


Diffusion  Coefficient,  cm  /sec 


Diffusion  Coefficient,  cm  /sec 


Figure  5 


p 

Plot  of  the  diffusion  coefficient  (cm  /sec)  vs  the  density  gradient 
scale  length  (cm)  for  various  values  of  the  neutral  density,  (cmj 
in?a  singly  ionized  and  neutral  oxygen  plasma,  n  =  10°  xm  ,  L  = 
10  cm,  B  =  0.35  gauss,  T  /T.  =  1,  and  p.  =  4.62  x  10^  cm.  the 
corresponding  classical  diffusion  coefficient  from  interspecies  cou 
lomb  collisions  is  1.68  x  104  cnr/sec. 


j * 
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Diffusion  Coefficient,  cm  /sec 


Density  Gradient  Scale  Length,  cm 


2 

Figure  6.  Plot  of  the  diffusion  coefficient  (cm  /sec)  vs  the  density  gradient 
scale  length  tan)  for  various  values  of  the  neutral  density  (cm" J) 
in,a  singly  ionized  and  neutral  oxygen  plasma,  n  *  107  xm  ,  L  = 
10  cm,  B  =  0.35  gauss,  T  /T.  =  1,  and  p.  =  4.62  x  10  cm.  the 
corresponding  classical  drffdii on « coefficient  from  interspecies 
coulomb  collisions  is  1.68  x  105  cm  /sec. 


Diffusion  Coefficient,  cm  /sec 


Figure 


Density  Gradient  Scale  Length,  cm 


2 

7.  Plot  of  the  diffusion  coefficient  (cm  /sec)  vs.  the  density 
gradient  scale  length  (cm)  for  various  values  of  the  neutral 
density  (cttr3)  in  a  singly  ionized  and  neutral  oxygen  plasma. 
0  =  10®  cm'3  l2  =  107  cm,  B  *  0.35  gauss,  WTj  =  1 ,  and 
*  4.62  x  102  cm.  The  corresponding  classical  diffusion  coef¬ 
ficient  from  interspecies  coulomb  collisions  is  1.68  x  10 
cn^/sec. 
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Diffusion  Rate,  sec 


Diffusion  Rate,  sec 


Diffusion  Rate*  sec 


Figure  10 


.  Plot  of  the  diffusion  rate  (sec  *)  (the  diffusion  coefficient 
divided  by  L  )  vs  the  density  gradient _fcale  length  (cm)  for 
various  values  of  the  neutral  density  (on  )  in  ajingly  ionized 
and  neutral  oxygen  plasma,  n  •  10  cm-  ,  L  »  10'  cm,  B  »  0.35 
gauss,  T./T.  ■  1,  and  p.  -  4.62  x  104  cmr  The  corresponding 


Diffusion  Rate,  sec 


•  ■  mimr  n"’”»lMlfT^^iiitf|i|ifir^ 


Density  Gradient  Scale  Length,  cm 


Figure  11.  Plot  of  the  diffusion  rate  (sec"  )  (the  diffusion  coefficient 
divided  by  l  )  vs  the  density  gradient  fcale  length  (cm)  for 
various  values  of  the  neutral  density  (on  )  in  a  iSingly  ionized 
and  neutral  oxygen  plasma,  n  =  10  cm  »,  L  *  10  ct,  B  =  0.35 
gauss,  T  /T.  =  1,  and  p.  «  4.62  x  10  cm.  The  corresponding 


E 


Diffusion  Rate,  sec 


Figure  12.  Plot  of  the  diffusion  rate  (sec-1)  (the  diffusion  coefficient 
divided  by  Lp2)  vs.  the  density  gradient  scale  length  (cm)  for 
various  values  of  the  neutral  density  (cm-3)  in  a  singly  ion¬ 
ized  and  neutral  oxygen  plasma,  n  =  108  cm"3.  Lz  =  107  cm, 

B  =  0.35  gauss,  Te/T.j  =  1,  and  p-j  =  4.62  x  102  cm. 
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103  104  105  io6 

Density  Gradient  Scale  Length,  cm 


"  2 
Figure  13.  Plot  of  the  diffusion  coefficient  (cm  /sec)  vs  the  density 

gradient  scale  length  (cm)  for  various  values  of  the  neutral 
density  (cmlj)  in  a  singly  ionized  and  neutral  oxygen  plasma, 
n  *  l(r  cm  ,  L  *  10°  cm,  B  *  .035  gauss,  p,  ■  4.62  x  Kr  cm. 
The  corresponding  classical  diffusion  coefficient  from  inter¬ 
species  coulomb  collisions  is  1.68  x  10^  crrr/sec. 
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Diffusion  Coefficient,  cm  /sec 


m 


2 

Figure  14.  Plot  of  the  diffusion  coefficient  (cm  /sec)  vs  the  density 
gradient  scale  length  (cm)  for  various  values  of  the  neutral 
density  (cmZ'3)  in  a  singly  ionized  and  neutral  oxygen  plasma, 
n  =  10°  cm  ,  L  =  10°  cm,  B  =  .035  gauss,  p.  =  4.62  x  10J  cm. 
The  corresponding  classical  diffusion  coefficient  from  inter¬ 
species  coulomb  collisions  is  1.68  x  10°  cnr/sec. 
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Diffusion  Coefficient,  cm  /sec 


Diffusion  Coefficient,  cm  /sec 


Density  Gradient  Scale  Length,  cm 


Figure  16.  Plot  of  the  diffusion  coefficient  (cm^/sec)  vs.  the  density 
gradient  scale  length  (cm)  for  various  values  of  the  neutral 
density  (cm-3)  in  a  singly  ionized  and  neutral  oxygen  plasma 
n  =  10'  cm-3,  lz  =  108  cm,  B  =  0.035  gauss,  p-j  =  4.62  x  103 
cm.  The  corresponding  classical  diffusion  coefficient  from 
interspecies  coulomb  collisions  is  1.68  x  107  cm^/sec. 


Diffusion  Rate,  sec 


Diffusion  Rate,  sec 


Figure  20.  Plot  of  the  diffusion  rate  (sec-^)  vs.  the  density  gradient 
scale  length  (cm)  for  various  values  of  the  neutral  density 
(cm"3)  in  a  singly  Ionized  and  neutral  oxygen  plasma,  n  * 
107  cm-3,  U  *  108  cm,  B  =  0.35  gauss,  ■  4.62  x  103  cm. 


Density  Gradient  Scale  Length,  cm 


Figure  22.  Plot  of  the  diffusion  coefficient  (cmz/sec)  vs.  the  density 
gradient  scale  length  (cm)  for  various  values  of  the  neutral 
density  (cm-3)  in  a  singly  ionized  and  neutral  barium  plasma 
n  =  10'  cm-3,  L,  =  2.5  x  106  cm,  B  =  0.35  gauss,  Te/T-j  =  1, 
and  pj  *1.1  x  TO3  cm.  There  is  no  anomalous  diffusion  at 
nn  =  10'°  cm-3  and  above.  The  classical  coulomb  diffusion 
coefficient  is  1.68  x  105  cmz/sec. 
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Figure  23. 


Density  Gradient  Scale  Length,  cm 


Plot  of  the  diffusion  rate  (sec-1) (the  diffusion  coefficient 
divided  by  L  i  vs  the  density  gradient  scale  length  (cm)  for 
various  values  of  the  neutral  density  (cnfg)  in,a  singly  , 

ionized  and  neutral  barium  plasma,  n  =  10°  cm  ,  L  =  2.5  x  10 
B  =  0.35  gauss,  T  /T.  =  1,  and  p.  =  1.1  x  10J  cm.  There  is  no 
_ jirr...:!.  1  *  -  -  ,  ^10l  5  -3  .  . 


anomalous  diffusion  at  n_  =  101 


and  above. 


Diffusion  Rate,  sec 


Density  Gradient  Scale  Length,  cm 


Plot  of  the  diffusion  rate  (sec-1)  (the  diffusion  coefficient 
divided  by  Li2)  vs.  the  density  gradient  scale  length  (cm)  for 
various  values  of  the  neutral  density  (cm--3)  in  a  singly  ion¬ 
ized  and  neutral  barium  plasma,  n  =  10'  cm"3,  Lz  =  2.5  x  10b 
cm,  B  =  0.35  gauss,  Te/Ti  =  1,  and  pf  1.1  x  103  cm.  There  is 
no  anomalous  diffusion  at  n„  =  10'°  cm"3  and  above. 


Table  1.  Modal  parameters  associated  with  diffusion  for  n  =  106  cm’3  and 
neutral  density  =  10  cm  .  The  diffusion  coefficients  and  rates 
(DIFRTM)  for  this  table  are  plotted  in  Figures  5  and  10.  The 
frequency,  kz  and  k  values  of  the  mode  with  largest  diffusion 
contribution  for  given  ambient  parameters  are  given  as  FREQM,  KZM, 
and  KYM,  respectively.  The  corresponding  modal  growth  rate  (y)  is 
(DIFRTM)  x  (L^/p^),  with  p.  =  4.62  x  10^  cm. 


L 

1.000000E+03 
1.258925E+03 
1.584893E+03 
1. 995262E+03 
2. 511886E+03 
3. 162278E+03 
3. 981072E+03 
5.011872E+03 
6.309573E+03 
7.943282E+03 
1. 000000E+04 
1. 258925E+04 
1. 584893E+04 
1.995262E+04 
2.511886E+04 
3. 162278E+04 
3.981072E+04 
5.011872E+04 
6.309573E+04 
7.943282E+04 
1.000u00E+05 
1.258925E+05 
1.584893E+05 
1.995262E+05 
2. 511886E+05 
3. 162278E+05 
3.981072E+05 
5.011872E+05 
6. 309573E+05 
7.943282E+05 
1.000000E+06 


DIFRTM 

3.814634E+00 
2. 369956E+00 
1.464790E+00 
9. 099748E-01 
5.609202E-01 
3.433089E-01 
2.090178E-01 
1.270536E-01 
7.607813E-02 
4. 530903E-02 
2.692350E-02 
1.609415E-02 
9.401440E-03 
5. 384977E-03 
3.037790E-03 
1.729756E-03 
9. 774986E-04 
5.460831E-04 
2.817955E-04 
1. 141932E-04 
4. 312406E-05 
1.250024E-05 
0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


FREQM 

-9. 769753E+00 

-9.496501E+00 

-5 . 978324E+00 

-5.479202E+00 

-4. 269953E+00 

-3.682328E+01 

-2.627965E+00 

-2.014682E-00 

-1.820443E+00 

-1. 218509E+00 

-9.000448E-01 

-8.273793E-01 

-5.893355E-01 

-4. 381297E-01 

-2.825324E-01 

-2.886011E-01 

-1.919144E-01 

-1.401367E-01 

-1.352709E-01 

-1. 1597 13E-01 

-8.447109E-02 

-7.567366E-02 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


KZM 

1.066358E-05 
1.066358E-05 
8.033977E-06 
8.033977E-06 
6.052827E-06 
6. 052827E-06 
4.337902E-06 
4. 1 181 16E-06 
3.585815E-06 
2.483167E-06 
2.411565E-06 
2. 149995E-06 
1 . 591634E-06 
1. 224249E-06 
8.441201E-07 
1. 129290E-06 
8. 117573E-07 
6. 283185E-07 
6.283185E-07 
6.283185E-07 
6.283185E-07 
6.283185E-07 
0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


KYM 

1. 382213E-03 
1. 338603E-03 
1. 328175E-03 
1. 350158E-03 
1. 139464E-03 
1.2012 12E-03 
1 . 02612  IE-03 
1. 121778E-03 
9.699496E-04 
8. 386703E-04 
9.623929E-04 
8.422823E-04 
7. 371620E-04 
6.451611E-04 
5.646424E-04 
6. 968087E-04 
6. 172797E-04 
5.468276E-04 
4.844164E-04 
4.291284E-04 
5.695165E-04 
5. 10667 1 E-04 
0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


Table  2. 


Modal  parameters  associated  with  diffusion  for  n  =  106  cm-3  and 

?nrcoTi!.wnsitu  =  cm  *  diffusion  coefficients  and  rates 
(DIFRTM)  for  this  table  are  plotted  in  Figures  5  and  10.  The 
frequency,  k2  and  k  values  of  the  mode  with  largest  diffusion 
contribution  for  given  ambient  parameters  are  given  as  FREQM,  KZM, 

/nrr-£LM?  refPect i vel y .  The  corresponding  modal  growth  rate  (y)  is 
(DIFRTM)  x  (Lj/pj),  with  p.  =  4.62  x  102  cm. 


L 

1.000000E+03 

1.258925E+03 

1.584893E+03 

1.995262E+03 

2.511886E+03 

3. 162278E+03 

3.981072E+03 

5.011872E+03 

6.309573E+03 

7.943282E+03 

1.000000E+04 

1.258925E+04 

1.584893E+04 

1.995262E+04 

2.511886E+04 

3. 162278E+04 

3.981072E+04 

5.011872E+04 

6.309573E+04 

7. 943282 E+04 

1.000000E+05 

1.258925E+05 

1 . 584893E+05 

1.995262E+05 

2.511886E+05 

3. 162278E+05 

3.981072E+05 

5.011872E+05 

6.309573E+05 

7.943282E+05 

1 . 000000E+06 


DIFRTM 

3.812435E+00 

2.368332E+00 

1.463419E+00 

9.089246E-01 

5.601502E-01 

3.427019E-01 

2.085587E-01 

1.266694E-01 

7.581483E-02 

4. 510241E-02 

2.674409E-02 

1.597287E-02 

9. 309413E-03 

5.317152E-03 

2.982904E-03 

1.689643E-03 

9.461926E-04 

5.252151E-04 

2.691069E-04 

1.078633E-04 

3.577331E-05 

8.427287E-06 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


FREQM 

-9.768440E+00 
-9.495257E+00 
-5. 976889E+00 
-5.477855E+00 
-4.268416E+00 
-3.680882E+00 
-2.626207E+00 
-2.012980E+00 
-1.818750E+00 
-1.216663E+00 
-8.982056E-01 
-8.255790E-01 
-5.874334E-01 
-4. 361742E-01 
-2.804537E-01 
-2.866735E-01 
-1.898572E-01 
-1.588993E-01 
-1.334874E-01 
-1. 145967E-01 
-8.277380E-02 
-7.437796E-02 
0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


KZM 

1.066408E-05 

1.066408E-05 

8.034318E-06 

8.034318E-06 

6.053055E-06 

6.053055E-06 

4. 337902E-06 

4. 118116E-06 

3.585815E-06 

2.483167E-06 

2.411565E-06 

2. 149995E-06 

1.591634E-06 

1.224249E-06 

8.441201E-07 

1. 129290E-06 

8.117573E-07 

7.010836E-07 

6.283185E-07 

6.283185E-07 

6.283185E-07 

6.283185E-07 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


KYM 

1.382213E-03 

1. 338603E-03 

1.328175E-03 

1. 350158E-03 

1. 139464E-03 

1.201212E-03 

1.026121E-03 

1. 121778E-03 

9.699496E-04 

8.386703E-04 

9.623929E-04 

8. 422823 E- 04 

7.371620E-04 

6.451611E-04 

5.646424E-04 

6.968087E-04 

6. 172797E-04 

5.468276E-04 

4.844164E-04 

4.291284E-04 

5.695165E-04 

5. 106671E-04 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 
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Table  3.  Modal  parameters  associated  with  diffusion  for  n  =  10  cm"  and 
neutral  density  =  108  cm"  .  The  diffusion  coefficients  and  rates 
(DIFRTM)  for  this  table  are  plotted  in  Figures  5  and  10.  The 
frequency,  k  and  k  values  of  the  mode  with  largest  diffusion 
contribution  ror  given  ambient  parameters  are  given  as  FREQM,  KZM, 
and  KYM,  respectively.  The  corresponding  modal  growth  rate  (y)  is 
(DIFRTM)  x  (L  /p. ) ,  with  p.  =  4.62  x  10^  cm. 


L 

1 . 000000E+03 
1.258925E+03 
1.584893E+03 
1.995262E+03 
2.511886E+03 
3. 162278E+03 
3. 981072E+03 
5.011872E+03 
6.309573E+03 
7. 943282E+03 
1.000000E+04 
1.258925E+04 
1.584893E+04 
1.995262E+04 
2.511886E+04 
3. 162278E+04 
3.981072E+04 
5.011872E+04 
6.309573E+04 
7. 943282E+04 
1 . 000000E+04 
1.258925E+05 
1 . 584893E+05 
1.995262E+05 
2. 511886E+05 
3. 162278E+05 
3. 981072E+05 
5.011872E+05 
6. 309573E+05 
7.943282E+05 
1.000000E+06 


DIFRTM 

3.790454E+00 
2.352100E+00 
1.449715E+00 
8. 984281E+01 
5.524566E-01 
3. 366383E-01 
2.039740E-01 
1.228337E-01 
7. 318907E-02 
4.304561E-02 
2.495868E-02 
1.476939E-02 
8.400296E-03 
4.651611E-03 
2.504937E-03 
1. 308483E-03 
6.767454E-04 
3. 3i2950E-04 
1.525947E-04 
4.981477E-05 
4.219431E-06 
0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


FREQM 

-9. 755299E+00 
-9.482805E+00 
-5. 962516E+00 
-5.464354E+00 
-4.253012E+00 
-3. 666386E+00 
-2.608569E+00 
-1. 995874E+00 
-1.801739E+00 
-1. 198066E+00 
-8.796328E-01 
-8.074161E-01 
-5.681944E-01 
-4. 163627E-01 
-3. 166900E-01 
-2. 124721E-01 
-2.041036E-01 
-1.390763E-01 
-1. 157783E-01 
-1.012316E-01 
-8. 144010E-02 
0. 

0. 

0. 

0. 

0. 

'  0. 

0. 

0. 

0. 

0. 


KZM 

1 . 06691 IE-05 
1 . 06691 IE-05 
8. 037728E-06 
8.037728E-06 
6.055340E-06 
6.055340E-06 
4.337902E-06 
4. 1 181 16E-06 
3.585815E-06 
2.483167E-06 
2.411565E-06 
2. 149995E-06 
1.591634E-06 
1.224249E-06 
9.784000E-07 
7. 144664E-07 
9.226763E-07 
7.010836E-07 
6.283185E-07 
6.283185E-07 
6.283185E-07 
0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


KYM 

1.382213E-03 
1. 338603E-03 
1. 328175E-03 
1.350158E-03 
1. 169434E-03 
1.2012 12E-03 
1. 026121E-03 
1. 121778E-03 
9.699496E-04 
8. 386703E-04 
9.623929E-04 
8.422823E-04 
7.371620E-04 
6.451611E- 04 
5.646424E-04 
4.941727E-04 
6. 172797E-04 
5.468276E-04 
4.844164E-04 
4. 291284E-04 
3.801506E-04 
0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 
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Table  4.  Modal  parameters  associated  with  diffusion  for  n  =  10  cm  and 
neutral  density  109  cm  .  The  diffusion  coefficients  and  rates 
(DIFRTM)  for  this  table  are  plotted  in  Figures  5  and  10.  The 
frequency,  k  and  k  values  of  the  mode  with  largest  diffusion 
contributionzfor  given  ambient  parameters  are  given  as  FREQM,  KZM, 
and  KYM,  respectively.  The  corresponding  modal  growth  rate  (y)  is 
(DIFRTM)  x  (L/p.),  with  pi  =  4.62  x  10  cm. 


L 

DIFRTM 

FREQM 

KZM 

KYM 

1 . 000000E+03 

3.571184E+00 

-9.622375E+00 

1.071929E-05 

1 . 3822 13E-03 

1.258925E+03  . 

2. 190344E+00 

-9. 356805E+00 

1.071929E-05 

1.338603E-03 

1.584893E+03 

1.313225E+00 

-5.816218E+00 

8.071748E-06 

1.328175E-03 

1.995262E+03 

7.940079E-01 

-5. 326672E+00 

8.071748E-06 

1.350158E-03 

2.511886E+03 

4.761895E-01 

-4.095235E+00 

6.078115E-06 

1. 139464E-03 

3. 162278E+03 

2.766258E-01 

-3.517393E+00 

6.078115E-06 

1 . 2012 12E-03 

3.981072E+03 

1.588103E-01 

-2.425808E+00 

4. 337902 E-06 

1.026121E-03 

5.011872E+03 

8.711108E-02 

-2. 145190E+00 

3.760364E-06 

8. 765516E-04 

6.309573E+03 

4.768608E-02 

-1.623686E+00 

3. 58581&E-06 

9.699496E-04 

7.943282E+03 

2.396089E-02 

-1. 379444E+00 

3. 122317E-06 

8.386703E-04 

1 . 000000E+03 

1. 148859E-02 

-8.767534E-01 

2.205375E-06 

7.251592E-04 

1.258925E+04 

4.663466E-03 

-7.220698E-01 

1.958659E-06 

6.270114E-04 

1.584893E+04 

1.599977E-03 

-4. 650747E-01 

1.468000E-06 

5.421476E-04 

1.995262E+04 

2.850033E-04 

-3.005671E-01 

1. 144637E-06 

4.687699E-04 

2.511885E+04 

0. 

0. 

0. 

0. 

3.162278E+04 

0. 

0. 

0. 

0. 

3.981072E+04 

0. 

0. 

0. 

0. 

5.011872E+04 

0. 

0. 

0. 

0. 

6.309573E+04 

0. 

0. 

0. 

0. 

7.943282E+04 

0. 

0. 

0. 

0. 

1 . 000000E+05 

0. 

0. 

0. 

0. 

1.258925E+05 

0. 

0. 

0. 

0. 

1.584893E+05 

0. 

0. 

0. 

0. 

1.995262E+05 

0. 

0. 

0. 

0. 

2.511886E+05 

0. 

0. 

0. 

0. 

3. 162278E+05 

0. 

0. 

0. 

0. 

3.981072E+05 

0. 

0. 

0. 

0. 

5.011872E+05 

0. 

0. 

0. 

0. 

6.309573E+05 

0. 

0. 

0. 

0. 

7. 943282 E+05 

0. 

0. 

0. 

0. 

1 . 000000E+06 

0. 

0. 

0. 

0. 

47 


6  -3 

Table  5.  Modal  parameters  associated  with  diffusion  for  n  =  10  cm  and 

neutral  density  =  1010  cm  .  The  diffusion  coefficients  and  rates 
(DIFRTM)  for  this  table  are  plotted  in  Figures  5  and  10.  The 
frequency,  k  and  k  values  of  the  mode  with  largest  diffusion 
contributionzfor  given  ambient  parameters  are  given  as  FREQM,  KZM, 
and  KYM,  respectively.  The  corresponding  modal  growth  rate  (y)  is 
(DIFRTM)  x  (L^/p^),  with  p..  =  4.62  x  10z  cm. 

L  DIFRTM  FREQM  KZM  KYM 


1.000000E+03 
1.258925E+03 
1. 584893E+03 
1. 995262E+03 
2.511886E+03 
3. 162278E+03 
3.981072E+03 
5.011872E+03 
6.309573E+03 
7.943282E+03 
1.000000E+04 
1.258925E+04 
1.5848 93 E+04 
1. 995262E+04 
2.511886E+04 
3. 162278E+04 
3. 98 1072 E+04 
5.011872E+04 
6.309573E+04 
7. 94 3282 E+04 
1.000000E+05 
1.258925E+05 
1.584893E+05 
1. 995262E+05 
2.511886E+05 
3. 162278E+05 
3. 981072E+05 
5.011872E+05 
6.309573E+05 
7. 943282 E+05 
1.000000E+06 


1.702275E+00 
7.903521E-01 
2.961617E-02 
7. 100461E-02 
0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


-1.037543E+01 
-7. 753612E+00 
-5.481703E+00 
-3. 514090E+00 
0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


1. 121000E-05 
8.403553E-06 
6.299708E-06 
4.622389E-06 
0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


1 . 175676E-03 
9.452586E-04 
7.600003E-04 
6. 110503E-04 
0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 
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where  o>'  is  the  mode  oscillatory  frequency,  n  is  the  ambient  ionized  density 

o  2 

outside  of  the  striation  vg  =  c|Ex|/B,  b..  =  (k^p..)  ,  and  IQ  and  are  the 

modified  Bessel  functions  of  zero  and  first  order.  (See  Eq.  (B-14)  and  the 
preceding  argument  in  Appendix  B  for  details.)  This  effect  is  due  to  convec¬ 
tion  of  the  mode  out  from  a  central  region,  i.e.,  the  drift-dissipative  eigen¬ 
mode  is  treated  non-locally  in  the  x-direction,  the  direction  of  both  the 
density  gradient  and  the  background  electric  field.  The  expression  for  the 
diffusion  rate  is  assumed  to  be  given  by  Eq.  (7)  with 


D  *  Tf/Te  -  !„<!>, >]3/2 


but  the  growth  rate,  which  scales  inversely  as  the  gradient  scale  length  at 
the  location  around  which  the  mode  is  centered,  is  less  by  a  factor  of  order 
(Vi*)av/ve  even  the  additional  decrement  of  Eq.  (12).  (See  Eq.  B-12) 
in  Appendix  B  and  the  proceeding  argument  for  details.) 

Hence  diffusion  and  growth  are  inhibited  whenever  ( v-j*)av/ve  <<:  with 
growth  inhibited  at  least  as  rapidly  as  (vi*)av/ve»  and  with  diffusion  at 
least  as  rapidly  as  [(v  *)  /v  ]2. 

1  a  V  “ 
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6.  CHARACTERIZATION  OF  THE  DIFFUSION  COEFFICIENT 
FROM  DRIFT-DISSIPATIVE  MODES 
IN  COMPUTER  SIMULATIONS 


Values  for  D*.  the  anomalous  diffusion  coefficient  in  the  absence  of  a 
background  electric  field  parallel  to  the  ionized  density  gradient,  can  be 
obtained  from  the  graphs  of  Figures  ( 3) - (7 )  and  (13)-(16). 

From  Eqs.  (B-12),  (B-13)  and  (B-14)  with  b.  «  1  and 

u-  =  6 

Lx  vea 

the  suggested  form  for  the  instability  growth  rate  when 
|ve|  >  yon^y 
is 


Y  * 


Is  JL 

L  i  a 


(Vf) 


av 


/  r 

£ 

^  \ A  \f  ~1 

.25  -  /T  f  - 

0\  ve 

"  /  v.  '  vi 

"-L  e 

Hence,  a  necessary  condition  of  instability  is  then 


r4± 

\  n  j  vt. 


<  0.18 


Now  def i  ne 


03) 


a0  = 


-  c(£  •  Vn 1 
BlVnl 


a 


_  Tflv n| 


1  m^n 


In  view  of  (13),  as  a  lower  bound  on  the  anomalous  diffusion  coefficient  D, 
we  have 


D  =  D*  ,  oig  <  otj 


D  =  0 


a0  >  al 


04a) 

04b) 
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7.  OBSERVATIONAL  ASPECTS  OF  DRIFT-DISSIPATIVE  MODES 

7.1  ESTIMATION  OF  THE  VARIANCE  DUE  TO  DRIFT-DISSIPATIVE  MODES 

The  variance  is  related  to  an  integral  over  the  power  spectral  density. 
The  contribution  to  this  integral  from  drift-dissipative  modes  is  likely  to 
be  almost  directly  proportional  to  a  second  integral  for  the  root  mean  square 
of  the  density. gradient  due  to  drift-dissipative  modes.  Since  the  second 
integral  can  be  estimated  by  physical  arguments,  a  plausible  estimate  is 
possible  for  the  variance,  even  though  the  power  spectral  density  is  not 
known  in  detail. 

We  employ  the  definitions^  ' 

VvV  =  ^7<C|/dx  /dye  x 

for  the  2-D  power  spectral  density  and 

CO  oo 

aND2a//VVky>  dkxdky 

-OO  -00 

for  the  variance. 

First,  represent  the  drift-dissipative  mode  as  a  single  mode  which 
varies  harmonically  in  the  y-direction.  We  have 

N(x,y)  *  nj(ky)  e  r  .  (16) 

Define 


f  -  /dx/ 


x  ^  ik..x+ikyy 


dy  e 

0  0 

After  some  algebra,  we  have 

1 2  _  1  _  /|  >\  1 2 


N(x,y) 


4K„ 


fl  -  v  lnT(kP1  6ky ,ky' 


52 


and 


Vkx*V  =  ifi2/(47r2xY> 

This  results  in 

2  |n,(k;)|2 

0nd  =  ~j^r 

For  purposes  of  estimates  on  using  Eq.  (5),  we  have 


lnl  (k\PI  ~  k  L 


x  x 


where 


L  =  gradient  scale  length  for  region  where  mode  is  located. 


and,  as  noted  in  Section  5, 


Kx  "  (piLjJ 


Then, 


'0 


1 


ND 


2  ;  2  KX 


=  n 


2  pi 


3/2 


Kx  Lx 


L/X 


If  we  identify  X  with  an  outer-scale  length  LQ  (possibly  a  rod  radius, 
assuming  striations  to  be  rods  of  enhanced  density  in  the  plane  xB),  a 
plausible  limitation  for  Lx  is 


4p,  <  Lx  <  L„ 


Hence, 


2  p 


n  2  3/2 

i  2  n0  pi 


3/2 


~TTZ  <  °ND 


Lo 


4  2 

As  examples,  for  Lq  =  4  x  10  cm,  p.  =  4  x  10  cm: 
"o2  ,0'3  <  aN02  <  "(/  5  x  ,0'3 
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6  2 

and  for  Lg  =  2  x  10  cm,  =  4  x  10  cm: 
nQ  3x10  <  cr^p  <  ng 


7.2  MODIFICATION  FOR  DRIFT  WAVE  REPRESENTATION  AS  A  SERIES  OF  ky  MODES 
We  now  let 


N(x,y)  =  ^Tn^kp  e 
k' 

y 

Stabilization  is  expected  when 
VN  ~  nQ/L^ 


-ikp-Kxlxl+^(ky) 


We  estimate 


VN  = 


-ikp-Kx|x|+i<J>(kp 


k; 


2\i 


If  there  is  no  correlation  between  the  various  (J>(kp,  then  at  |x|  =  0 


2  „  2U 


VH£  |ni(k;)|  kx 


X' 

y 


If  K  does  not  vary  significantly  over  k',  as  seems  reasonable  for  drift- 
x  y 

dissipative  modes,  then 


E  i"i(kv 


"O 

2  2 
Lx  Kx 


In  parallel  with  the  single-mode  calculation. 


f  =EV  "!<kP  e 


-id)  (k-)  2K 


Kx2  +  kx2  v; 
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and 


_  2  1  Y"*  |  |2  _  1  0 

aN  *nLini(yi  “  kt rr ;  2 

x  k-  x  Kx  Li_ 

This  is  the  same  expression  as  for  the  single-mode  case,  so  that  the  same 
numerical  estimates  apply. 


7.3  ESTIMATE  OF  THE  MAXIMUM  PSD 

It  is  reasonable  to  assume  that  the  maximum  in  the  PSD  due  to  drift 

waves  is  peaked  at  kx  m  and  ky  m  corresponding  to  the  modal  parameters  for 

the  maximum  linear  growth  rate,  and  that  the  variance  can  be  calculated  by 

assuming  that  the  maximum  is  roughly  constant  over  two  regions  of  range 

Ak  ^  2k  ,  Ak  &  k  with  negligible  contribution  beyond, 
x  x,m  y  y,m 

Then, 


^N^x.m’^y.m^  ^x,m  ^y,m  ~  °N 
or  using  k^m  *  kx  =  (pjLJ"*,  ky>ro  =  .4/pr  results  in 

-4/Pi  ]  ^ -6  nQ2  p.3/X  . 


(17) 


(14) 

This  result  can  be  compared  against  the  two-dimensional  fornr  of 

the  isotropic  power  spectral  density  evaluated  at  (k  ,k  )  for  a  spectrum 

^2  ^  y 

which  in  one-dimension  varies  as  k  : 


4>(kx,ky)  = 


4tt^  An2  r(3/2)  L 
4tt2(L‘2  +  kx2  +  ky 


-1 


(18) 


(Here  and  below,  we  have  used  the  results  of  Ref.  (14)  for  <j>,  multiplied  by 

2  1 

an  additional  factor  (4ir  )  in  order  to  preserve 


yj(kx, ky)  dkx  dky  =  A?  .) 
— 1  2 

For  p^/L  ^  .16  and  An  =  nQ  / 4, 

<0(kx,ky)  *  .22  nQ2  p.3/L 


(19) 
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(20) 


2  2 

For  p./L  «  .16  and  An  =  nQ  /4, 

d>(kx,ky)  ~  .6  nQ2  p.3/L 

By  way  of  comparison,  the  contribution  from  a  cylindrical  rod  of  radius  L 
2  2  i 

■in  the  limit  (k  +  k  r  L  »  1,  taking  the  average  over  rapidly  varying 
x  y 

phases,  is  about  0.72  of  (18)  in  the  same  limit 


7.4  OBSERVABILITY  OF  THE  DRIFT  WAVE  PSD  THROUGH  IN-SITU  MEASUREMENT 

Comparison  of  (17)  with  (19)  or  (20)  indicates  the  possibility  of 

-2 

detection  of  drift  waves  against  the  2-D  equivalent  of  a  1-D  k  background. 


One  has 


n02  <  Vkx,m,ky,m) 

4^2  ~  ^=2Tkx,m’ky,m] 


2.5  n. 


where  4>D  denotes  the  theoretically  estimated  PSD  contribution  from  drift- 
dissipative  modes  and  denotes  the  contribution  from  the  2-D  counter¬ 
part  of  a  1-D  k-2  PSD. 

Comparison  (21)  is  appropriate  to  the  detectability  of  drift  waves 
through  in-situ  measurement.  To  see  this,  assume  a  probe  moving  in  the  x^- 
direction  and  no  variation  in  the  x3-direction.  The  measured  Fourier  com¬ 
ponent  is 


nl  (kl  ,(P 


/n0 


-iklxl 


(where  we  have  averaged  over  the  X2~direction) ,  and  the  spectral  density 
is 


Wh.iP  a  fni^ki,o^2 

where  the  subscript  1  refers  to  the  x^ coordinate  axes.  If  the  x^  axis 
makes  an  angle  0  with  the  x-axis  of  our  previous  analysis,  then 

♦m(kl  0>  =  Vkl  cos0’kl  sin0) 


Hence,  a  peak  in  the  2-D  PSD  could  be  observable  provided  k  ~  k.  cose, 

mx  1 

k  ~  k,  sin0.  There  are  three  caveats: 
my  l 
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1.  For  kmy  »  kmx,  the  range  of  0  with  appreciable  drift  mode  con¬ 
tributions  will  be  small,  if  in  addition  the  various  (localized) 
regions  with  drift  waves  are  not  randomly  oriented. 

2.  The  quantity  4>^(kl,0)  or  4>[sj(kxrn»kym)  is  an  average.  A  single  rod 

of  radius  L  can  (within  a  factor  of  order  unity)  account  for  a  1-D  PSD 

_2 

of  form  k  with  outer-scale  length  L.  If  the  drift  waves  are  found 
only  on  the  edges  of  a  small  number  of  large  rods,  one  could  have  an 
observation  pattern  in  which  most  measurements  showed  no  contribution 
with  an  occasional  very  large  contribution. 

3.  Even  if  >k  )/(()  _7  <  1,  the  anomalous  diffusion  coefficient 

could  still  be  much  larger  than  the  coulomb  diffusion  coefficient. 

7.5  OBSERVABILITY  OF  THE  DRIFT  WAVE  PSD  THROUGH  SCINTILLATION  MEASUREMENT 
We  now  consider  the  contribution  of  the  theoretically  derived  drift- 
dissipative  waves  to  scintillation,  especially  as  a  diagnostic  of  drift  wave 
behavior  for  the  correlation  functions  of  the  log  amplitude  and  the  phase 
fluctuations  in  the  Rytov  approximation. 

Paraphrasing  Ref.  (2),  one  has  respectively: 


,  r”  iik,  r\  . 

x(*,)  x(v^  =  ^J  sr  J  **  s'" 

J  _rr>  I  1 


2k  (x2s  +  x2g'  x2^ 


4>2(x2,k1 )  exp(kl]p) 
(23) 


_______  2  f ”  dk1  fXs 

.(X.,)  <f>s(xl+pT  =  k  J  —  J  .  dz 


cos 


l_2k  ^x2s  +  x2g  '  X2^J 


4>2(x2,k1)  exp(ik1p) 
(24) 


where 

k  =  rf  wave  number 

x2  =  coordinate  along  propagation  line  of  sight 
x-|  =  coordinate  perpendicular  to  line  of  sight 
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x2s  =  tur^u^ent  layer  thickness 

X£g  =  distance  from  bottom  of  turbulent  layer  to  receiver 

The  turbulent  layer  has  no  variation  in  the  third  direction  (parallel  to 
the  magnetic  field). 

In  Eqs.  (23)  and  (24), 

I  An .  (k-,  ,k~  =  0)  j  2  ? 

^2 (x2 * )  =  - X — X - =  ^  4>(k-j  »0)  a 

where  An.  is  the  deviation  of  the  index  of  refraction  from  its  mean,  $  is 
1  2 

the  2-D  density  PSD,  and  a  is  the  proportionality  constant  between  density 

variation  and  index  of  refraction  variation.  For  w  »  u>  (with  to  the  elec- 

2?  P  P 

tron  plasma  frequency),  a  =  4rre  /(mew  ). 

For  drift  waves  with  >  kxm,  a  necessary  condition  for  scintillation 
contribution  from  a  given  in-situ  £,  typified  by  (kxm,k^m)  is  e2  •  t  =  0, 
where  is  the  unit  vector  in  the  X2~direction.  Therefore,  drift  waves  con¬ 
tribute  to  scintillation  only  if  the  magnitude  of  the  angle  0  of  rf-wave  pro¬ 
pagation  with  respect  to  the  x-direction  (i.e.,  the  direction  of  the  back¬ 
ground  density  gradient  projected  perpendicular  to  the  ambient  magnetic  field) 

*  *  _i 

is  less  than  0  where  9  =  tan  (^xm,k  ).  The  range  in  magnitude  for  con¬ 
tributing  k1  is  k  .  Hence,  (with  ex>  the  unit  vector  in  the  x-direction)  from 
Eqs.  (17)  and  (25), 


.  _  ’6  "o'  pi  „_2„  2 

2  ^2*^1 '  ~  ^ 


k, 

pi  '  pi 


cos-1 (k~-e  / I kp I }  <  tan'’(k-m/k,„)  • 


'  <L  x  1  ^ 1  '  xm  ym' 

'bp ( ^2 » k-j )  =  0  otherwise.  (2/ 

Denote  drift  wave  contributions  by  the  subscript  D  and  power  law  con- 
tributions  from  a  2-D  k  PSD  by  the  subscript  P.  For 

x2s  +  x2g/x2  <  k  <  ^x2s  +  x2gH-08)/pi2  (2? 

from  Eqs.  (23)  and  (24),  we  have 

X02(x1  )/Xp2(x1 )  >  <i>D2(x])/(|)p2(x1) 
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Given  Eqs.  (23),  (26),  and  (28), 


XD2(Xl)  =  2.38  x  10 2(-^)2  n-0-"f  --2s 

u  1  Wr/  vr  x 


From  Eq.  (24a)  of  Ref.  (14)  for  X£g  k.  X2$» 


Xp2(Xl)  -  3.62  x  101  Z-4-)2  a7  ^2-^. 

y  1  \mc  f  r  X 

Therefore,  the  condition  for  drift  wave  dominance  is 


k  >  k* .  -X  *:£  _L 

-  K  6.57  _  2  2 


nQ  Pi 

2  * 
The  maximum  value  of  Xg  (x^ ) ,  with  Eq.  (29)  satisfied,  occurs  with  k  =  k 

Then, 


Xg  (x-,)  =  1.0 


x  104(4) 
'roc  / 


2\2  V  6  X2s 

7  ZHT~  Pi  - T~ 

:)  (An2)2  XX2g 


c  o 

_  As  an  example,  using  the  data  of  Ref.  (15)  (e.g.,  nQ  =  6  x  10  cm"  , 

An^  =  3.6  x  10^  cm  ,  pi  =  3.35  x  102  cm,  X£S  =  2  x  10^  cm,  X ^  =  3  x  10^  cm, 
x  =  2.  x  106  cm),  one  has: 

k  -  .41  cm’^ 

Xg2(Xl)  =  4.6  x  10"5 

Hence,  there  can  be  rather  stringent  conditions  on  the  observability  of  drift 
waves  by  scintillation. 

7.6  ESTIMATE  OF  DECAY  OF  THE  PRIMARY  (LARGE  WAVELENGTH)  VARIANCE  IN  TERMS 
OF  THE  DRIFT  WAVE  VARIANCE 

It  is  clear  that  Eq.  (8)  applies  for  a  diffusion  coefficient  Dx  due  to 
drift-dissipative  modes.  On  using  Eqs.  (5)  and  (7),  Eq.  (8)  becomes 


&/"*  «•-*  /<*> 


k/  n,2  dV 
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where  the  right-hand  side  of  the  equation  is  due  to  drift  wave  contributions 
On  neglecting  variations  in  kx,  we  have 

2 

2 

TT  =  "2  <Y>  aND 

2  ? 
where  is  the  larger  wavelength,  non-drift  wave  variance  and  o^D  is  the 

drift  wave  variance.  Hence,  for  a  time  scale  At  for  decay  of  a we  have 


At  ~ 


L±{cm) 


aND2  8.6  x  103 


Here,  is  the  scale  length  for  density  variation  in  the  region  where  drift 
waves  are  found. 

We  caution  that  Eq.  (30)  is  equivalent  to  a  time  derivative  evaluated 

2 

at  t  =  0.  Variations  in  and  <y>  with  increasing  time  will  affect  the 
time  for  a.,2  to  decrease  a  given  non-infini tessimal  amount. 


8.  DISCUSSION  OF  EXPERIMENTAL  DATA  WITH  RESPECT  TO 
DRIFT  WAVES  AND  DIFFUSION 


8.1  8UARO 

The  BUARO  barium  plasma  injection  experiment  conducted  by  Group  J -10  of 
the  Los  Alamos  Scientific  Laboratory  in  June  1976  provides  some  data  which  may 
be  directly  relevant  to  drift  mode  diffusion  of  striations.  Figure  6  of  the 
Quick-Look  Report, (our  Figure  25),  ranges  from  about  250  km  in  altitude 
(bottom  of  photograph)  to  400  km  in  altitude  (top  of  photograph).  The  stria¬ 
tions  are  oriented  along  the  magnetic  field  lines  which  make  an  angle  of  about 
45 ii  with  the  vertical.  It  would  appear  that  the  cloud  is  optically  thin  espe¬ 
cially  at  its  higher  altitudes.  The  photograph,  taken  at  810  seconds  after 
the  injection,  indicates  more  diffuse  structure  at  higher  altitudes  in  that 
the  gaps  between  striations  at  lower  altitudes  appear  to  fill  in  at  higher  al¬ 
titudes.  More  quantitative  measurements  bear  this  out.^17^ 

We  consider  the  possibility  that  the  diffusion  at  higher  altitudes  is 
due  to  drift-dissipative  modes  which  are  not  as  strongly  unstable  at  lower  al¬ 
titudes.  Using  a  CIRA  Model  2  atmosphere  at  4  AM^8^  the  neutral  number  densi- 

9-3  7-3 

ties  at  250  km  and  400  km  are  approximately  10  cm  and  2  x  10  cm  .  The 

barium  ion  column  density  has  been  estimated  at  about  8  x  1011  cm-2,^7^  lead- 

4  -3 

ing  to  a  representative  barium  ion  density  of  about  B  x  10  cm  .  The  column 
density  of  ambient  ionization  is  roughly  the  same.  We  choose  ion  and  electron 
temperatures  of  0.1  and  0.2  eV,  respectively.  These  values  are  somewhat  arbi¬ 
trary  and  will  be  discussed  below.  Since  the  barium  photoelectrons  are  emit¬ 
ted  with  an  energy  of  order  0.6  eV  and  since  the  barium  ions  can  have  an 
energy  of  order  1.5  eV  it  is  reasonable  to  choose  temperatures  above  the  ambi¬ 
ent  neutral  temperature  of  0.07  eV. 

Results  for  the  diffusion  coefficient  as  a  function  of  L  for  barium  are 

5  6 

shown  in  Fig.  (22)  for  n  =  10  and  n  the  neutral  density,  equal  to  10 

in  n 

through  10  with  a  scale  length  parallel  to  the  magnetic  field  of  100  km. 
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Density  Gradient  Scale  Length,  cm 


2 

Figure  26.  Plot  of  diffusion  coefficient  (iji  cm  Vsec)  vs  L  (in  cm)  for?sin 
gle  ionized  barium  ions,  n  =  ICr  cm"'5,  T.  =  O.F  eV,  L  =  10/  cm 
B  =  0.35  gauss.  z 


The  suppression  of  diffusion  for  n  =  10'^  and  L,  >  0.1  km  when  com- 
q  n 

pared  to  nn  <  10  is  qualitatively  in  keeping  with  the  possibility  that  drift 
dissipative  modes  are  suppressed  at  lower  altitudes  for  the  BUARO  release 

9 

(although  the  actual  densities  at  250  km  and  400  km  in  the  release  are  10  and 

7-3  9-3 

2  x  10  cm  respectively).  For  n_  <10  cm  ,  the  diffusion  coefficients,  0 

**  -j  o  6 

,  as  calculated  are  somewhat  below  10'  ctrr/sec  and  the  region  with  D  i  >  10 

2  9-3 

cm  /sec  extends  up  to  L  ,  =  0.3  km.  Thus  for  n  <  10  cm  as  contrasted  to  n 

10  _3  **■  n  " 
=10  cm  ,  over  a  time  of  1000  seconds  the  calculations  indicate  significant 

diffusion  for  scale  lengths  in  the  0.1  -  0.5  km  range. 

We  note  however  that  diffusion  if  it  affects  the  barium  cloud  features 
is  more  likely  to  be  in  the  1.0  -  5.0  range. 

The  comparison  is  of  course  only  semi-quantitative.  An  increase  in  the 
barium  ion  on  neutral  collision  frequency  of  a  factor  of  three  would,  roughly 
speaking,  lower  the  neutral  density  appropriate  to  each  contour  by  a  factor  of 
three.  The  values  of  the  diffusion  coefficient  would  be  increased  by  larger 
temperatures  and  possibly  by  more  accurate  eigenmode  calculations,  especially 
of  kx.  Since  the  ion  dynamics  of  the  modes  is  primarily  in  the  perpendicular 
direction,  any  noise  with  wavelength  primarily  parallel  to  the  magnetic  field 
could  act  to  increase  the  electron  collision  frequency  without  similarly  in¬ 
creasing  the  ion  collision  frequency.  With  the  background  parameters  we  use 
this  would  encourage  growth. 

Further  differentiation  between  growth  at  higher  and  lower  altitudes 
would  occur  if  there  is  significant  barium  ion  temperature  variation  as  a 
function  of  altitude.  Barium  ions  can  lose  energy  through  collisions  with 
atomic  oxygen  ions  (time  scale  of  35  seconds  with  ionized  oxygen  density  of 
105  cm-3).  The  oxygen  ions  can  then  engage  in  resonance  charge  exchange^173) 
with  oxygen  atoms  (time  scale  of  about  5  sec  at  300  km,  and  about  50  sec  at 
400  km).  Hence  assuming  an  e- folding  time  for  barium  temperature  decay  of  85 
sec  at  400  km  altitude  it  would  appear  that  barium  ion  temperatures  of  greater 
than  0.2  eV  can  exist  there  for  a  time  of  order  several  hundred  seconds.  The 
temperature  loss  process  would  be  slower  if  any  appreciable  fraction  of  am¬ 
bient  oxygen  ions  recombined  with  explosion  byproducts. 

8.2  STRESS 

Data  from  the  rocket  probe  experiments  for  STRESS  exhibits  high  fre¬ 
quency  behavior,  specifically  a  peak  in  ( Ar»RMS/n)  as  a  function  of  frequency. 
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(Here  AnR^<.  is  the  root  mean  square  ionized  density  irregularity  amplitude.) 
Together  with  other  rocket  parameters,  data  for  the  peak  can  be  summarized  as: 


Rocket 

Velocity 

Event 

Frequency 

Perpendicul ar  to  B 

Parallel  to  B 

Esther 

=  900  sec"1 

0.8  km/ sec 

1.04  km/sec 

Fern 

-  1200  sec"1 

0.26  km/sec 

1.34  km/sec 

The  peak  is  observed  prior  to  entering  the  barium  cloud  (Case  A)  (our  Fig. 

27a,  Fig.  6.1  of  Ref.  19),  and  while  in  the  barium  cloud  but  not  in  the  stri¬ 
ated  region  (Case  B)  (our  Fig.  27b,  Fig.  6.3  of  Ref.  19).  For  the  barium 

cloud  striated  region  there  is  a  flat  region  of  (AnRM<-/n)  from  about  10  Hz  up 

to  the  frequency  designating  the  peaks,  then  a  sharp  drop  above  the  peak  fre¬ 
quency  (Case  C)  (our  Fig.  27c,  Fig.  6.7  of  Ref.  19).  The  value  of  (AnRMS/n) 
appears  appreciably  higher  outside  of  the  barium  cloud  than  in  the  striated 
region  of  the  barium  cloud. 

It  is  likely  that  the  peak  in  (AnRMS/n)  represents  density  fluctuations 
parallel  to  the  magnetic  field.  For  fluctuations  parallel  to  B  the  wavelength 
is  given  by 

rocket  velocity  parallel  to  B 
frequency 

For  fluctuations  perpendicular  to  B  the  wavelength  is  given  by 

rocket  velocity  perpendicular  to  B 

frequency  ' 


Figure  27.  Irregularity  Amplitude  vs  Frequency  for  the  Esther  Release  of 
Project  Stress  as  Measured  by  Probe  Rocket  at  Release  Plus 


If  the  density  fluctuation  is  parallel  to  B  the  wavelength  is  about  1.1  meter 

for  both  Esther  and  Fern.  If  the  density  fluctuation  is  perpendicular  to  B  it 

is  about  0.9  meter  for  Esther  and  0.2  meter  for  Fern.  Hence  if  the  wave  has 

the  same  wavelength  in  both  cases  it  is  parallel  to  B. 

Further,  the  value  of  (AnRM<Vn)  appears  to  decrease  within  the  striated 

region  compared  to  either  of  the  other  two  regions.  From  Figure  (5.3)  of  Ref.  9, 

2  3 

it  would  not  be  surprising  if  An  in  the  range  of  10  -10  Hz  were  roughly  the 
same  for  Figures  27a  and  27c.  For  modes  dependent  on  density  gradients,  one 
would  expect  (AnR^<./n),  if  anything,  to  increase  the  striated  regions,  since  the 
density  gradients  are  larger  there.  In  addition,  for  the  altitudes  (~  170  km) 
and  density  (ng  >  10  cm”  )  of  the  Project  Stress  barium  clouds,  theoretically 
we  do  not  expect  short  wavelength  drift  wave  phenomena  to  be  significant. 

These  considerations  lead  us  to  believe  that  the  high  frequency  dis¬ 
turbances  are  not  related  to  drift  mode  diffusion  perpendicular  to  the  magnet¬ 
ic  field.  The  difference  between  Cases  B  and  C  is  suggestive  of  nonlinear 

(19) 

plasma  behavior.  The  existence  of  density  scale  lengths  of  order  15  meters' 
is  consistent  with  the  absence  of  diffusive  instabilities  on  the  sides  of 
striations  as  expected  theoretically  at  these  altitudes  and  densities. 

8.3  THREE  METER  BACKSCATTER 

Backscatter  by  3  meter  irregularities  moving  vertically  with  vertical  k 
vector  is  frequently  observed  in  Equatorial  Spread  F.^2(^  One  possible  ex¬ 
planation  for  this  phenomenon  is  that  the  3  meter  disturbance  is  either  an 
unstable  drift-dissipative  mode  or  is  driven  by  somewhat  larger  wavelength 
unstable  drift-dissipative  modes.  The  proposed  geometry  is  shown  in  Figure 

28.  Along  the  sides  of  the  depleted  region  there  may  be  a  background  electric 

(21 ) 

field,  Ex,  in  the  direction  of  the.  density  gradient.  ' 

In  this  subsection  we  treat  the  connection  between  drift-dissipative 
modes  at  3  meter  wavelengths  and  the  power  spectral  density  at  longer  wave¬ 
lengths,  which  would  be  expected  to  dominate  the  contributions  to  anomalous 
diffusion.  A  power  law  fall-off  for  the  ID  PSD  of  the  form  k”  ,  v  =  7  would 
be  required  to  link  both  the  theoretically  calculated  PSD  at  ky  =  0.4/p^  and 
the  observed  PSD  at  ky  =  2.1  x  10"2  cm. 
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gure  28.  Schematic  geometry  for  d 
backscatter.  (Wavel°ngth 


We  first  present  a  detailed  treatment  of  backscatter  contributions  from 
power  law  power  spectral  density,  then  determine  the  power  law  variation,  and 
then  present  some  discussion  (based  on  other  theoretical  calculations)  of 
factors  which  would  work  in  the  direction  of  a  rapid  fall-off  of  the  PSD  with 
increasing  k. 

The  formula  for  the  P„ynting  flux  from  backscatter  is: 

2  2 

<S(r)>  s  w  (h)  Re  JJ  n(r')l  A0(r')|  2  sin2xOO  B£(r',r") 

x{exp  1 k[ ( r  -  r'\  -  |r  -  r"|  +  |R  -  r'\  -  |R  -  r" | ]} 

(221 

All  quantities  are  as  defined  in  Tatarskii.  For  our  purposes  it  is  suffi¬ 
cient  to  specify 

Be(r',r")  =  ex(r")  e^r") 


where  e,  is  the  dielectric  constant.  Further  for  id  z  2 id  (with  <d  the  elec- 
1  pe  pe 

tron  plasma  frequency) 


(31) 


with  i  the  sum  over  electrons.  We  note 

y^5(r>  -  r.)  d3r  =  n(r')  d3r 
i 

with  n  the  electron  density 

Using  (31)  we  may  write: 

2 

e1(r')e1(r")  *  IX)5js*r’  "  ri'5(r"  '  ; C r~  -  r()S(r“  -  rj) 

\  »  «/  i  rfi 

The  first  double  summation  is  the  contribution  from  incoherent  backscatter; 
the  second  double  summation  is  the  contribution  from  coherent  backscatter. 
Hence  the  ratio  of  coherent  backscatter  to  incoherent  backscatter  is 
I^kyi^k)  where: 

I2(k)  =  ff ’  ri^r"  ‘  ri )  exp{ik[|r  -  r'|  -  |r  -  r"| 

V  V  i  i'^i 
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|R  -  r"|  -  jR  -  r" i  ] }  d3K  d3r"  , 


Ijfk)  =  ffY.  &{r'  -  r.)6(r"  -  r^)  exp{1k[|r  -  r'\  -  |r  -  r" 


V  V  i 


+  |R  -  r'\  -  |R  -  r"  |  ]}  d3r'  d3r"  . 

23  24 

to  within  about  a  factor  of  2  accuracy:  ’ 
Ij(k)  -  n  /d3r' 


where  V  is  the  backscatter  volume.  For  backscatter  at  ky  =  k^ : 

1  C  f  ~  9  ik  x-ik.x2/r  ik  z-(ik./r)z2 

l2  =  ^2  J  J  ln(kx  -  2Vkz>l  e  e  dkxdxdkzdz 

(32) 


n(lc)  *  y*n 


-ller  A  v 
e  d3r  . 


[r  is  the  distance  from  the  scattering  volume  to  the  receiver  (sender).] 

2 

For  a  power  law,  power  spectral  density  with  variance  <6n>  and  ID 

-v  14 

power  spectral  density  varying  as  k  ,  we  have 

ln(kx»ky»kz)l2  =  “  “  ~  r~~v/2+l  (33) 

(L  +  kX  +  ky  +  Ve  ) 


A  *  <6n>^X2r(v/2  ±  l.).(2w)3_  #  (34) 

4*  Lv_1  r(3/2)  r[(v  -  D/2] 

(z  is  the  direction  of  the  magnetic  field  with  anisotropy  oarameteri zed  by  e.) 

-  2 

The  normalization  for  Ini  follows  from  the  constraint 


i(£)|  2  d3k  =  (2 


x3  /  2,+ 
»)  / n  (r 


(r)  d2r  =  V(2tt)3  <6n>2 


We  define  z*  to  be  the  spatial  extent  of  the  scattering  volume  in  the 
z-direction.  Then  three  convenient  limits  obtain  from  Eqs.  (32)  and  (33): 
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Case  1 


„  i  [*h 

**  ~  qi  •  V  — 


z*  «  1  ; 


Case  2 


W  =  2n  A 


i  rd/2)  rr(v  ±  ii/r 


7  O  (v/2)+(l/2)  2  r[(v/2)  +  1. 

[1/r  +  (2k!)  ] 


1  fih 

kje  ’  V  irr  z*  1  ; 


_ l _ i  r(i/2) 

.  „  -  (v/2)+(l/2 )  2  r[( 

[1/L2  +  (2k!)2] 


rf(v  ±  l ) /2 1  rSE 

v/2)  +  1] 


Case  3 


2*  >> 


=  7  ’  2  (v/2)+l  2  *  (37) 

[1/L2  +  (2k!)2] 

Woodman  and  Basu15  have  examined  the  compatibility  of  a  one-dimensional 
power  spectral  density  with  v  =  2  and  the  observed  ratio: 

I9(k)  ?  1 

log  tStt  *  40  db  ,  k  =  2k.  =  2.09  x  10'*  cm"1  .  (38) 

1 

In  the  interest  of  clarity,  we  repeat  the  evaluation  of  log[I7(k)/I, (k)]  using 

6  c  A  4  -3 

their  experimentally  determined  values:  L  =  2  x  10  cm,  <Sn  =  6  x  10  cm  ,  n 
=  6  x  108  cm  3,  Pj  =  3.35  x  102  cm. 

If  one  associates  the  3  meter  behavior  with  drift-dissipative  modes,  by 

_2 

our  numerical  calculations  the  most  unstable  kz  at  ky  =  1.04  x  10  cm,  n  = 
105  cm-3,  n  =  108  cm"3  is  k  ~  10"5  cm"*.  One  then  has 


Then  from  "Case  3"  with  v  =  2  we  obtain 


(~)  = 


84  db  . 
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This  compares  with  the  previously  estimated33  value  of  82  db.  Thus  is  appears 
likely  that  the  v  =  2  spectrum  does  not  persist  down  to  backscatter 
wavelengths. 

Alternately  one  can  ask  for  the  value  of  v  consistent  with  the  theoret¬ 
ical  drift-dissipative  wave  estimate  of  Eq.  (17)  at  ku  =  0.4/p.  and  with  the 
.  “  _?  Li 

experimental  observation  of  Eq.  (39)  at  =  2kj  =2.1x10  cm  . 

Experimentally,  using  (37)  and  (38): 


W 


A 

e 


_ I_ 

(2k1) 


V 


(40) 


Let  us  assume  that  4>N ( 'cXm*kyni) *  as  given  (17)>  and  ttie  three  meter  back¬ 
scatter  behavior  are  linked  by  the  same  power  law. 

On  using  (17),  (33)  and  (34)  we  obtain: 


^N^xm’^ynP 


<6n> 


4it  L 


& 


II 


1 


(-2  2  2  \( v/2 )+ (1/2 ) 

L  L  +  r  +  r 
1  xm  ym / 


0.6  n2P3L 
o  1 


(41) 


On  using  (41)  to  eliminate  <6n>2/Lv'3  in  A  in  (40)  and  taking  L"2  «  k(L  + 
j  xm 

k  ,  we  obtain  for  k  =  k  =  0.4/p.,  v  *  7.3  and  for  k  =  0,  k  -  0.4/p., 
ym  xm  ym  i  xm  ’  ym  i  ’ 

v  »  6.4. 


Hence  the  theoretically  calculated  PSD  at  k  m,k  is  not  necessarily 

xm  ym  J 

inconsistent  with  the  3  meter  observations  but  a  very  rapid  fall-off  with 
increasing  k  is  required.  Such  behavior  could  occur  if  smaller  values  of  n 
and/or  L  are  required  for  instability  to  break  out  at  short  wavelengths. 

Our  numerical  calculations  with  drift-dissipative  modes  tend  to  bear 
this  out  (although  lower  hybrid  modes  have  not  been  examined  in  detail).  The 
following  observations  are  available  for  drift-dissipative  modes  (DDM): 

1.  At  3  meter  wavelengths  there  is  no  instability  from  DDM  for  any  n 
>  104  cm"3,  106  <  nn  <  1010  cm'3,  l_x  >  103  cm  for  Tg  =  T .  =  0.1  eV,  B  =  0.35 
g,  Lz  =  107  cm. 

2.  For  wavelengths  in  the  range  6-12  meters,  there  can  be  instability 

for  n  <  105  cm"3,  n  <  10^  cm"3. 

~  n  *  4^3  8-3 

3.  For  the  particular  case  of  n  =  10  cm  ,  n  *  10  cm  :  at  6  meters 

3  4  n 

there  is  instability  for  10  <  Lx<  1.25  x  10  ,  at  12  meters  there  is  insta¬ 

bility  for  103  <  L,  <  4.0  x  104  cm. 

5  _3  8-3 

4.  For  the  case  of  n  =  10  cm  "  ,  np  =  10  cm"  :  at  6  meters  there  is 
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3  3 

stability,  at  12  meters  there  is  instability  for  10  <  L^<  5.0  x  10  cm. 

5.  On  the  other  hand,  for  ODM  including  all  wavelengths  with  kyLj_ >  1, 
from  Figures  3-7,  one  has  for: 

n  =  104,  nn  =  108,  instability  for  103  <  Lj_  <  4.0  x  104 

n  =  108,  nn  =  108,  instability  for  103  <  L±  <  8.0  x  104 

n  =  108,  np  =  108,  instability  for  IQ3  <  <  108 

On  the  basis  of  these  calculations.it  appears  possible  that  stable  DDM 

disturbances  at  3  meters  could  be  driven  nonlinearly  by  unstable  DDM  at  6 

meters,  and  that  it  would  not  be  unreasonable  for  this  to  happen  preferen- 

4-3  5 

tially  at  lower  background  plasma  densities  (of  order  10  cm  rather  than  10 
-3  6  3 

cm  or  10  cm--1).  Hence  the  very  low  PSD  at  3  meters  would  result  from  the 

combination  of  factors: 

1.  The  mode  at  3  meters  is  stable. 

2.  The  ambient  density  of  the  localized  region  in  which  the  nonlinear 
driving  takes  place  may  be  relatively  low  compared  to  the  average  ambient 
density  (104  cm-3  vs  108-108). 

3.  If  the  ambient  density  is  not  relatively  low  compared  to  the  aver¬ 
age  ambient  density,  then  the  volume  occupied  by  regions  of  the  then  requisite 
high  density  gradients  could  be  very  small  or  non-existent. 

As  an  additional  aside  we  note  that  the  presence  of  a  range  of  excited 
ky  at  the  same  altitude  and  time  should  result  in  Frequency  Spread  F. 
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9.  CONCLUDING  REMARKS 


1.  It  appears  possible  that  plasma  turbulence  involving  modes  other 
than  the  gradient  drift,  E  x  B  mode^  can  contribute  significantly  to 
striation  dynamics  at  high  altitudes.  The  BUARO  barium  plasma  release  at  450 
km  altitude  has  a  diffusion  pattern  which  is  qualitatively  similar  to  that 
expected  from  drift-dissipative  mode  turbulence.  Three  meter  backscatter  in 
Equatorial  Spread  F  can  be  explained  as  due  to  drift-dissipative  modes  [e.g., 
stable  modes  being  driven  by  nearly  (in  t  space)  unstable  modes].  In  both 
cases  other  explanations  are  possible;  however,  the  explanations  involving 
drift  dissipative  modes  have  been  followed  through  in  some  detail  and  within 
current  observational  and  theoretical  uncertainties  appear  to  be  consistent 
with  the  phenomena. 

2.  There  are  several  aspects  to  the  role  of  drift  dissipative  modes  in 
diffusion  and  the  power  density  spectrum. 

At  striation  tips  the  outbreak  of  drift  modes  can  effectively  transfer 
fluctuations  over  a  range  in  length  scale  typically  from  0.5  km  to  0.05  km. 
Drift-dissipative  waves  should  furnish  small  scale  diffusion  of  large  scale 
density  gradients  and  inhibit  bifurcation,  provided  the  maximum  length  scale 
of  the  modes  producing  diffusion  is  much  smaller  than  the  bifurcation  wave¬ 
length,  e.g.,  Eq.  (10)  is  satisfied,  and  provided  the  resulting  diffusion  is 
sufficiently  large  to  satisfy  a  criterion  of  the  form  of  Eq.  (9a)  [or  perhaps 
Eq.  (9b)].  Since  drift-dissipative  wave  diffusion  itself  is  a  sensitive  func¬ 
tion  of  L^,  bifurcation,  in  practice,  will  be  limited  by  the  smaller  of  the 
value  of  satisfying  condition  (11)  or  the  value  of  (characteristically 
of  order  0.5  km)  for  which  D  becomes  sufficiently  large  to  satisfy  (9a). 
(This,  of  course,  is  based  on  the  assumption  v^n  »  t"1  where  \k  is  the  ion- 
neutral  collision  frequency  and  -t  is  the  time  scale  for  plasma  velocity 
change.)  The  suppression  of  bifurcation  does  not  necessarily  prevent  the 
outer  scale  length  of  a  striation  (=  s,Q)  from  subsequently  decreasing  from 
additional  plasma  convective  behavior  but  it  appears  likely  to  be  a  control¬ 
ling  influence  on  i  in  many  cases  of  interest.  The  value  of  0.5  km  for  t  is 

0  (251  0 

consistent  with  data  for  outer  scale  lengths  in  Equatorial  Spread  F. 
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The  disposition  of  the  drift  dissipative  modes  can  probably  occur 
through  nonlinear  processes  such  as  the  production  of  long  wavelength  vortex 
modes  as  well  as  through  linear  diffusion  processes. 

Either  direct  nonlinear  estimates^  or  wave-particle  diffusion  in 
drift-dissipative  turbulence^)  result  in  estimates  for  the  power  density 
spectrum 

2  p. 

~~n~  =  ~  » 

where  n..(k)  is  the  modal  amplitude  for  density  as  a  function  of  £.  This  spec- 
trum  is  flat  as  a  function  of  k+  therefore  looking  different  from  the  fre- 
quently  measured  k  *"  power  density  spectrum.  The  occurrence  of  the  spectrum 
at  striation  tips  which  occupy  a  relatively  small  fraction  of  the  perimeter  of 
a  striation  does  not  contradict  the  existence  of  drift-dissipative  modes  at 
striation  tips. 

Along  the  sides  of  striations  drift-dissipative  modes  can  act  to 
diffuse  plasma  and  presumably  alter  the  power  density  spectrum.  Except  at  its 
tip,  a  single  striation  can  be  taken  to  be  elongated  in  the  y-di recti  on  with 
density  n(x)  and  Ex(x)  related  by  the  condition  of  continuity  of  the  Pedersen 
current  in  the  x-direction.  This  is  valid  provided  the  ion-neutral  collision 
frequency  is  larger  than  the  reciprocal  of  the  time  scale  associated  with 
plasma  inertial  effects. 

It  is  shown  in  Appendix  B  and  Section  5  that  drift-dissipative  waves 
are  inhibited  whenever 

^  »  Pi[vti/(cEx/B)]  .  (42) 

(Here  ^  is  a  typical  scale  length  for  the  x-directed  density  gradient.) 
Therefore  drift-dissipative  wave  effects  are  most  likely  in  regions  across  the 
thin  (x)  dimension  of  a  striation  whereof  is  small.  If  the  inner  scale  length 
for  the  power  density  spectrum  corresponds  to  a  real  region  in  coordinate 
space,  drift-dissipative  wave  turbulence  will  tend  to  develop  in  this  region 
and  we  would  expect  3  meter  backscatter  to  originate  from  such  regions.  The 
observation  of  inner  scale  lengths  of  order  tens  of  meters  in  Equatorial 
Spread  F^25)  is  consistent  with  (42)  if  one  takes  *  3  meters  and 

Cvti/(cEx/B]  “  10* 

Drift-dissipative  waves  in  such  a  case  would  not  have  to  alter  a  k"^ 
power  density  spectrum  except  by  affecting  the  inner  scale  length  and  spatial 
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frequencies  above  the  spatial  frequency  corresponding  to  the  inner  scale 
length.  The  effect  on  diffusion  and  scintillation  decay  would  depend  strongly 
on  the  nonlinear  behavior  and  coupling  of  the  modes. 

3.  Calculation  of  drift-dissipative  growth  rates  under  model  HANE  con¬ 
ditions  as  exemplified  in  Figures  3-7  indicates  that  diffusion  peaks  at  n  < 
10  cm  for  the  parameters  chosen.  At  higher  ionized  densities  the  ion- ion 
collision  frequency  becomes  so  large  that  long  wavelength,  low  frequency  modes 
are  the  only  ones  possible.  Electron-ion  collisions  and  the  resulting 
ambipolar  diffusion  further  dampen  the  waves.  Increasing  neutral  density  acts 
monotonically  to  inhibit  the  waves. 

Results  for  a  model  of  a  "heaved1'  plasma  are  shown  in  Figures  13-20. 
Here,  too,  at  higher  altitudes,  lower  magnetic  field  values  and  larger  grad¬ 
ient  scale  lengths  parallel  to  the  ambient  magnetic  field,  turbulent  diffu¬ 
sion  can  be  significant. 

4.  On  the  basis  of  this  report  we  suggest  a  number  of  experimental  and 
observational  approaches  with  theoretical  interaction  for  examining  striation 
decay,  scintillation  behavior  and  the  occurrence  of  drift  modes: 

(a)  In-situ  measurement  is  likely  to  be  the  key  to  direct 
observation  of  drift-modes  at  the  relatively  long  wavelengths  (=  60 
meters,  for  an  0+  plasma  with  B  =  0.35  gauss)  expected  to  be  primarily 
responsible  for  anomalous  diffusion.  The  contribution  of  drift  modes  to 
scintillation  appears  to  be  weak  in  comparison  to  a  k  one-dimensional 
power  spectral  density. 

(b)  For  drift-dissipative  waves  (or  even  gradient-drift  waves)  the 
flux  due  to  turbulent  diffusion  can  be  approximated  as 


where  n^  and  ^  are  respectively  the  density  and  electric  field 
associated  with  a  plasma  mode,  and  the  brackets,  <>,  denote  a  suitable 
space  and  time  average.  Hence  in  situ  electric  field  measurements  and 
determination  of  (43),  which  depends  on  the  correlation  between  n^  and 
Ej,  would  be  desirable.  Expression  (43)  is  a  vector  expression;  ideally 
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one  would  want  to  know  its  direction  relative  to  the  background  (zero 
order)  plasma  density  gradient.  (If  E^  is  parallel  to  the  background 
density  gradient,  no  diffusion  results.) 

(c)  It  is  possible  that  barium  cloud  releases  in  the  200-300  km 
altitude  range,  for  which  some  drift  dissipative  wave  turbulence  is 
expected  (see  Figs.  21-24)  could  exhibit  drift-dissipative  wave 
phenomena.  Since  barium  cloud  density  gradients  can  be  determined  rather 
reliably  optically,  it  would  seem  that  one  could  obtain  a  good  idea  of  the 
direction  of  expression  (43)  with  respect  to  the  zero-order  plasma  density 
gradient. 

O 

(d)  Computer  simulations  of  the  type  performed  by  NRL,  using 
diffusion  coefficients  as  defined  in  Section  6  would  be  useful  in 
determining  and/or  verifying  observational  consequences  of 
formation  due  to  anomalous  diffusion.  (This  would  apply  to  structure  at 
length  scales  typically  greater  than  or  equal  to  the  drift  wave  length 
scales.)  Such  simulations  would  also  be  appropriate  for  assessing 
observational  differences  between  barium  cloud  striation  behavvior  in 
the  200-300  km  range  and  that  of  releases  at  lower  altitudes. 

(e)  Rocket  probe  measurements  of  density  alone  as  a  function  of 
wave  number  in  Equatorial  Spread  F  could  aim  at  establishing 
correlations  between  plasma  turbulence  and  the  decay  of  Spread  F  both 
along  the  elongated  sides  of  Spread  F  density  structures  and  at  their 
tips.  Additionally  if  drift  mode  turbulence  is  present  at  striation 
tips  we  would  expect  the  horizontal  widths  (or  outer  scale  lengths)  of 
ionization  structures  (similar  to  individual  striations  in  barium 
clouds)  to  have  a  minimum  size  of  order  0.5  km,  and  it  would  appear 
relatively  simple  to  determine  this  with  rocket  probes.  (Similar 
measurements  could  be  done  for  Mid-Latitude  Spread  F.) 

(f)  Backscatter  measurements  can  determine  the  power  spectral 
density  at  wavelengths  short  compared  to  the  expected  maximum  drift- 
dissipative  contribution  at  about  60  meters.  If  measurements  could  be 
made  at  in-situ  wavelengths  greater  than  three  meters,  some 
extrapolation  up  to  the  60  meter  wavelength  range  might  be  possible. 
The  three  meter  disturbances  if  they  are  drift  waves  are  not  as 
significant  for  diffusion  as  longer  wavelength  drift  modes.  Larger 
wavelength  drift  waves  are  limited  less  by  ion-ion  diffusion  and  hence 
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should  occur  earlier  in  time.  The  measurements  could  also  be  correlated 
with  scintillation  decay. 

A  range  of  vertical  wave  vectors,  k^,  at  the  same  altitude  and  time 
should  result  in  frequency  Spread  F. 

(g)  Usage  of  a  mid-latitude-shape  charge  release  aimed  vertically 
upwards  at  an  altitude  sufficiently  high  that  the  mean  free  path  is 
greater  than  the  scale  height  with  particle  velocities  of  order  5  km/sec 
would  serve  to  create  striations  under  conditions  more  likely  to  be 
characteristic  of  HANE  conditions,  with  ionized  density  greater  than  the 
ambient  at  altitudes  substantially  greater  than  200  km.  We  expect  drift- 
dissipative  turbulence  to  be  prominent  under  these  conditions  (see 
Figs.  13-20). 
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APPENDIX  A 


PHYSICAL  MECHANISM  FOR  THE  GROWTH  OF 
THE  ORIFT  DISSIPATIVE  MODE 

The  growth  of  the  drift  dissipative  mode  requires  that  the  wave  fre¬ 
quency  have  a  smaller  magnitude  than  the  diamagnetic  drift  frequency,  i.e., 

(i)  ^  *  (A-l) 

Provided  the  full  electron  and  ion  dynamics  permit  Eq.  (A-l)  to  be  satisfied, 
we  show  in  this  appendix,  using  simple  physical  arguments,  that  electron  col¬ 
lisions  destabilize  the  drift  dissipative  mode. 

The  electron  dynamics  of  the  drift  dissipative  mode  are  described  by 
the  following  set  of  fluid  equations 


0  =  -iunj  +  3(nvx^)/8x  +  3(nvz^)3z 

(A-2a) 

vxl  =  -icyi/B 

(A-2b) 

ik 

'll'.,.  <«n*l  -  "lTe> 
e  e 

(A-2c) 

with  first  order  variables  denoted  by  the  subscript  "one."  Equation  (A-2a) 
states  that  the  rate  of  change  of  the  perturbed  electron  density  per  unit  time 
is  balanced  by  the  flux  of  electrons  in  the  ambient  gradient  direction,  x,  and 
the  magnetic  field  direction,  z.  The  velocity  associated  with  the  x-directed 
flux  is  the  E  x  B  velocity  in  the  mode  (Eq.  (A-2b)).  The  velocity  in  the 
z-direction  is  a  consequence  of  the  balance  of  electron  collisions  with  the 
wave  electric  field  and  pressure  (Eq.  (A-2b)). 

Substituting  Eqs.  (A-2b)  and  (A-2c)  into  Eq.  (A-2a)  results  in  the  ex¬ 
pression  for  the  perturbed  electron  density  previously  calculated^) 

nl  .  q*l  /iu)DE  '  *\  (A-3) 

n  '  Te  \  ioi  -  «  ) 

(Note  that  in  this  paper  temperature  is  multiplied  by  Boltzmann's  constant.) 


By  using  Eq.  (A-3),  Eq.  (A-2c)  can  be  rewritten  in  the  following  way: 


ik. 


'  zl 


m  v  n 
e  e 


,  .  (»  -  kvvnr)(w  -  ifi) 

(qn*j)  - 


'l  7  7" 

0)  +  0 


(A-4) 


The  work  done  by  the  wave  is  then  equal  to 
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<4Ezl 


jzlEzl^  = 


_a!n 


to(« 


kyvDE} 


2  .  .2 

(0+6 


(A-5) 


Hence,  if  Eq.  (A-l)  is  satisfied,  work  is  done  by  the  plasma  on  the  mode  and 
the  drift  dissipative  mode  becomes  unstable.  A  qualitative  understanding  of 
the  instability  can  be  obtained  by  recognizing  that  if  Eq.  (A-l)  is  satisfied 
and  the  growth  rate  is  sufficiently  small  one  has: 

( nj/n |  >  q^/Tg  .  (A-6) 


Then  along  the  magnetic  field  lines  the  electron  flow  from  diffusion  is  larger 
than  the  electron  flow  from  mobility,  and  the  electrons  tend  to  flow  in  a  di¬ 
rection  opposite  to  the  electric  force  on  them.  In  this  case  the  electrons 
lose  energy  to  the  wave  and  the  wave  grows. 
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APPENDIX  B 


ELECTRIC  FIELD  LIMITATION  OF  DRIFT  DISSIPATIVE 
MODE  DIFFUSION 


The  purpose  of  this  calculation  is  to  evaluate  the  effect  of  an  electric 
field  Ex( x)  on  drift  dissipative  mode  growth  rates  and  diffusion.  We  assume 
that  the  background  density  gradient  is  in  the  x-direction  and  that  the  mode 
is  harmonic  in  y  and  z. 

The  motivation  for  this  study  is  the  observation  that  striations  perpen¬ 
dicular  to  the  magnetic  field  tend  to  be  elongated  in  the  direction  of  the  E  * 
B  drift.  Hence  around  most  of  the  perimeter  of  a  striation  the  electric  field 
should  be  parallel  (or  anti-parallel)  to  the  density  gradient.  [£x(x)  i  0.] 
[With  Ex(x)  =  0  previous  work  shows  that  drift  mode  growth  rates  are  likely  to 
be  large  enough  to  lead  to  diffusion  of  density  profiles  in  the  0. 1-1.0  km 
range. 

For  collisionless  ions  and  electrons  with  no  ambient  electric  field  the 

(26' 

drift  wave  equation  has  the  form'  ' 

2~ 

^4  +  Q ( x ,oj)  ?  «  0  .  (B-l ) 

d  L 

where 

Q  -  i.jl  ♦  [V  .*(*)]  Xoaj/Ta  .  (B-2) 


with 


t  =  i 


2r 

pLu) 


a  a 


«*(x)]  (X 

Oa 


-  X.  )/T 

la  a 


1  /  2 

Here  a  denotes  the  species,  v+=  =  (T  /m  )  ,  p  =  v.  /u  ,  w*  =  k.J-K/m.w  . , 

la  a  a  a  La  Ca  y  G  1  C  i 

0)e*  =  -  (1)*,  (i)j*  =  u*Tj/Te,  K  =  -  (dn/dx)/n,  and  L^  =  1/|K|  is  the  density 
scale  length  at  x  =  0,  the  point  around  which  the  mode  is  centered.  Further: 


with  b 

a 

^  » i 


Xj,a  =  W  eXP^ba)  Z kZVta>^  kZVta  * 


2  2  (27) 

=  ky  Pa  and  Z  the  plasma  dispersion  function.  '  We  note  for 


(B-3 ) 
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with  o=0,  Im  5  >  0;  o  =  1,  Im  5  =  0;  o  =  2,  Im  £  <  0. 

We  now  modify  Eq.  (B-l)  to  account  for  drift-dissipative  modes  in  the 
presence  of  the  electric  field,  Ex(x): 

(a)  In  the  limit,  k^v^/u^  <<  the  contribution  from 

ji  *  ue* 

is  equal  to: ^ 

q^n 

as  calculated  for  electrons  of  charge  -q.  The  corresponding  quantity  in  the 
limit  ve  »  a),  kzvte  is^10) 


with  6  =  kz2Te/meve  and  VDE  *  (+c/qB)TeK. 

(b)  The  summation  over  a  in  a  involves  finite  Larmor  radius  contribu¬ 
tions.  Since  we  expect  finite  electron  Larmor  radius  effects  to  be  small  we 
take 

A  -  A.  H  p2[*  +  «*(x)](Xoi  -  X1.)/T1  •  . 

(c)  Further  we  assume  that  the  sole  effect  of  the  electric  field  Ex(x) 
is  to  introduce  a  Doppler  shift  so  that  the  wave  frequency  o>  is  replaced  by  u> 
+  kyCEx/8  (=  O  wherever  it  occurs.  Defining  vg  =  c|Exl/B  we  neglect  differ¬ 
ences  in  the  gyrofrequency  and  gyroradius  of  order  (v  /v  .) (p./L  .)  and  correc- 
tions  to  the  ion  orbits  at  frequency  2u>ci  of  order  (v  /vt.)(p./L^)  with  re¬ 
spect  to  the  gyroradius. 


With  these  assumptions,  Eq.  (B-l)  can  be  approximated: 


i-i  +  0lx><01  =  0  (B-5) 

dx^  Ai 

with 

Q(x,«0  =  jl  +  «'  +  “f(x)]Xoi|/Ti  +  (ikyvDE  -  6)/(1«  -  6)Te  (B-6) 

To  illustrate  the  effect  of  the  electric  field  it  appears  sufficient  to 
consider  drift  modes  in  the  limit  s  »  kv^,  u  =  kv^.  Then  around  an  arbi¬ 
trarily  chosen  origin  one  may  expand 

*2 

Q(x)  =  Q(o)  +  xQ'(o)  +  f-  Q~(o)  .  (B-7 ) 

The  local  theory^0  ^  corresponds  to  setting  (d^$/dx^)  =  0  and  choosing 
the  origin  so  that  Q'(o)  =  0.  If  the  degree  of  localization, 

a  =  2(2)1/2  |Q~(o)f1/4  ,  (B-3) 

is  much  larger  than  the  wavelength  in  the  y  direction,  i.e.,  i  »  xy,  the 
local  theory  can  be  expected  to  be  approximately  satisfied.  Of  course,  the 
nonlocal  theory  to  be  discussed  subsequently  gives  an  even  more  accurate  de¬ 
scription  of  the  mode. 

Using  Eqs.  (B-3)  and  (B-6)  one  obtains 


The  condition  dQ/dx 


We  have: 


=  0  results  in 


d2q  =  _ 1 _  /  dV  b  A?\ 

dx2  TiAi“'  \  dx2  dx2/ 
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The  ion  momentum  equation  has  the  form: 


qv.  x  B 


o  =  — i - v.$.+  .a£_L_ 

m-c  it  m.  n 


7n 


from  which  to  lowest  order  in  v^/w  . : 


v  =  _  Il_  dn\_j_ 

IX  \  m.  m.n  dx/  2 

Vi  i  /  «ci 


In  the  limit  ve/“ce  <<  Vi/“C1-  the  condition  of  quasi-neutrality  results  in 


-2—  (n  v.  )  =  0 
ax  '  ix' 


from  which 


Hence 


and 


as  well  as 


'  +  *  .  .  u  cExo  no 
1  +  - « +  ky  -g—  fr 


(B-10) 
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Tlx  '  T^r["  (1  +  Ti/Te)“3)T  "  a  ky  ~T  ~  n  cTx J 


dx^  T^A^oo'L 


d  cd*  cE  „ 

(1  ♦  T-j/Tg)  V  9  ky  H®- 


"037(737)]  • 


(B-ll) 


For  Exq  =  0,  dQ/dx  *  0  is  equivalent  to  du^/dx  -  0  and  the  mode  is  centered 
around  an  extremum  of  l/n(dn/dx).  This  is  the  situation  conventionally  anal¬ 
yzed.  For  £xo  j  0,  from  dQ/dx  =  0  and  Q  =  0,  one  has: 


T. 

1 


m.u  • 
i  ci 


2 

1  d^n 
dx 


0 


The  first  term  in  the  equation  is  of  order  v^*(l/n) (dn/dx)  for  a  and  T./T 

order  unity.  (We  note  v..*  =  u^*/ky.) 

If  v  »  v4*  one  has 
e  i 


1  jin 

n  dx 


Vci 


of 


the  condition  for  mode  localization  is  not  met.  Hence  we  anticipate  the  mode 

2  2  1/2 

to  be  located  where  the  value  l/n(d  n/dx  )  '  is  of  order  the  average  value 
found  over  the  entire  region  of  non-zero  density  gradient.  If  is  the  typi¬ 
cal  scale  length  for  the  region  one  has: 


1_  d2n  _  B__ 
n  dx2  ~c*2 


with  e  of  order  unity  and: 


i  dn  ■  ~^i  S  1 

Vci  7  ave 


This  is  the  same  as: 


<*!>» 


-Ti 


m. o>  -ok 

1  ci 


an  average  value  for  the  ion-diamagnetic  drift  velocity  in  the  region  of  den¬ 
sity  gradient  one  has 


1  _  1_  ^vi^av  ^ 

rr "  se  a 


(B-12) 


Hence  the  mode  is  located  where  the  density  gradient  (L"1)  is  appreciably  less 

1 

[by  a  factor  of  (v.*).„/vj  than  the  mean  (=5^  )  and  its  growth  rate,  propor- 

^  ,  i  av  G 

tional  to  L  ,  is  correspondingly  decreased. 

For  v„  *  v.*,  we  expect  L,  =  . 

e  1  -*-(■)  o  1 

On  the  basis  of  previous  work'  '  instability  for  the  drift-dissipative 
mode  is  characterized  by 

cs 

rf  >  4  (B-13) 

J_  i 

1/2  2 
where  c$  is  the  ion  sound  speed  =  (Tg/m^)  and  v..  =  vin  +  0.3  i  ( kyPi )  , 

with  vin  the  ion-neutral  collision  frequency,  the  ion-ion  collision  fre¬ 
quency  and  the  deviation  of  Tg/Tj  from  unity  is  neglected.  For  vg  »  (v1-*)av 
a  necessary  condition  for  instability  becomes 
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s 
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4v 

TvTT 


i '  av 


a 

8 


or  equivalently: 


X  < 


The  effect  on  the  diffusion  coefficient  Dj_  can  be  estimated  through  the 
formula 


with  1/k^2  »  2/2  | [Q'^Co)]1^2]-1.  For  vg  &  v^*,  there  is  no  significant  modi¬ 
fication  in  the  previous  results  for  Q"(o)(^)  which  yields  kxo(^f 

For  v  »  v.*  one  has  -  given  that 
e  i 
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1  d2n 
n^  dx^ 
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and  using  Eqs.  (B-3),  (B-4),  (B-6)  with  Q  =  0,  (B-9a)  with  dQ/dx 
( B— 11 ) : 


Ti  -bi'13 

W 


tW  -  ■!»>,)]  « 


ve 

TWT 


i  av 


1 


p  .j<£ 


0,  and 


(B-14) 


We  now  discuss  the  effect  of  non-localness  of  the  mode. 

From  the  forms  of  Eqs.  (B-5),  (B-7)  and  (B-14)  it  is  apparent  that  the 
mode  has  a  nonlocal  mode  structure  about  x  =  0.  As  shown  in  Figure  29  an 
analogous  quantum  mechanical  picture,  that  the  mode  is  centered  a  the  top  of  a 
potential  hill,  can  also  be  made.  The  precise  solution  of  Eq.  (B-5)  can  be 
written  in  terms  of  Dp(z),  the  parabolic  cylinder  functions^28  ^  of  order  p, 
where 

°i?+7+p-^)  vz) 

For 


z  =  x/i. 

Q(o)  =  (7  +  p)  *  2 

q~(o)  =  - 1  r4 

and  a  boundary  condition  corresponding  to  the  outward  propagation  of  energy, 
the  appropriate  eigenmode  solution  to  Eq.  (B-5)  is 

♦  (x)  =  ♦q  Dp[exp(-iir/4)  |2Q"(o)|1/4  xj 

or  for  | x |  ♦  ® 
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♦  ( x)  =  <t>0  exp(-ipir/4)  |Q"(o)|p/4  xp  x  exp  i 
where  is  a  constant.  Eigenvalues  are  given  by  the  relation 


Q ( o )  =  -  i(2)1/2  (£  +  p)  |Q~(o)|1/2 


p=  0,  1,  2,  .  . 


Assuming  the  damping  rate,  y,  is  much  smaller  than  the  wave  frequency,  then 
from  Eqs.  (8-9). and  (8-14) 


Y  *  -  (1  +  2p)  71 


1  +  T. 
+  2d)  71  L - 1 


-  Il) 


-b.ru 

n  -  I  e  1  ,n  ,1/2  v 
e  o  /  o\  e 

77/7  /b'7  7i  “ 


(8-15 ) 


From  the  form  of  Eq.  (B-13)  it  is  evident  that  the  minimum  value  of  the  damp¬ 
ing  rate  occurs  for  p  =  0.  It  also  appears  that  electric  field  gradients  act 
like  magnetic  shear  to  stabilize  drift  modes  by  convective  damping. ^29^  Elec¬ 
trostatic  ion-cyclotron  , waves  propagating  in  the  ion  diamagnetic  drift  direc¬ 
tion  also  tend  to  be  damped  by  wave  convection  in  the  x-direction. ^ 

The  physical  mechanism  leading  to  the  damping  decrement  given  by  the 
expression,  Eq.  (8-15),  is  the  convection  of  wave  energy  away  from  x  =  0. 
Roughly  then  the  damping  decrement  must  scale  as  the  group  velocity  in  the 
x-direction  divided  by  the  scale  length  in  the  x-direction,  i: 


Then  using  the  replacement 


(9u  \  1 


i  k*(x) 


in  Eq.  (B-5)  results  in  an  expression  for  the  damping  rate  which  is  in  agree¬ 
ment  with  Eq.  (B-15)  to  within  a  constant  and  so  verifies  the  convective  ex¬ 
planation  for  the  damping. 


93 


APPENDIX  C 


l 


DERIVATION  Or  ‘10DAL  GROWTH  RATES  TAKING  INTO  ACCOUNT 
FINITE  ION  LARMOR  RADIUS  AND  ELECTRON-ION  COLLISIONAL  EFFECTS 


Drift  modes  with  frequencies  below  the  ion-cyclotron  frequency 
(i.e.,  w  <  £2.)  can  be  driven  unstable  in  an  inhomogeneous  plasma  if  the 
wave  frequency  is  smaller  than  the  diamagnetic  drift  frequency, 


e.g. , 


(0  <  (Dr 


.hA 

ai  Ln 


where 


(C-l) 


vi 


=  Ti/mi 


fii  = 


B0/mi 


(C-2) 


In  Eqs.  (C-l)  and  (C-2),  the  quantities  e,  c,  k  ,  and  Ln  are  the  magnitude 

of  the  electron  charge,  the  speed  of  light,  the  y-directed  wavenumber,  and 

the  density  gradient  scale  length,  respectively. 

In  a  collisionless  plasma,  the  mechanism  responsible  for  the  insta- 

f  51 

bility  is  inverse  electron  Landau  damping'  '.  Electron  diffusion  parallel 
to  the  magnetic  field  drives  the  instability  in  a  collisional  plasma.  In 
either  case,  it  is  apparent  from  Eq.  (C-l)  that  low-frequency  drift  insta¬ 
bilities  may  exist  in  plasmas  with 

kyP^  «  1  ,  Pi/Ln  <<  1  (C-3) 

where 

p,-  =  VS’  (C-4) 

is  the  ion  Larmor  radius. 

The  dispersion  relation  for  the  low-frequency  can  be  derived  using 
the  linearized  kinetic  equation,  Eq.  (IV. 100),  in  Ref.  5. 


(Vs-iW+  ikzVz-fis  gg)  fls(v)  = 


9*  3fn  (r,v)  n 

i  ■ —  k  • - — - +  v  — —  f .  ( r ,  v ) 

ms  9v  s  n0s  Os'  ’  ' 


1  s 


(C-5) 


where 
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MOW*  F**  BUMK-NOT  BUffl 


fls<M.t>  =  fls(t)  exp[i (k^y  +  k£z  -  wt)] 
$1 (r,t)  =  exp[i (kyy  +  kzz  -  a>t)] 
n1$(r,t)  =  nls  expCHk^  +  k?z  -  wt)] 


(C-6a) 

(C-6b) 

(C-6c) 

(C-6d) 

(C-6e) 


In  Eqs.  (C-5)  and  (C-6),  subscript  "0",  subscript  "1",  and  subscript  "s” 
denote  the  ambient  quantities,  perturbed  quantities,  and  arbitrary  species 
designations.  Also,  ^  is  the  perturbed  potential.  The  symbol  v$  repre¬ 
sents  the  effective  collision  frequency,  x,  for  species  s. 

A  straightforward  approach  to  the  solution  of  Eq.  (C-5)  is  based  on 
a  Fourier-Bessel  analysis.  We  write 


fls(v)  «  exp li  sin  (e)1y] exp(-ipe)ffs(vx,v11 

L  s  J  p=0  s 


vll)  (C-7) 


and  the  inverse  transformation 


fi5(vj.'vn>  -  ktT^l  s'  £  1 . 

L  s  j  o  s 

Here,  Jp  is  the  Bessel  function  of  the  first  kind  of  order  p.  One  finds 

'is  *  -  7-7!  <kzvz  +  P»s  -  “  +  «V'  Ih  *  kz»*Xl  +  rV*  -  “Dsl 
ms  vs  L  y  n  J 


(C-8) 


K 

j-3/2 

(2»v  2 

\-3/2 

\  s 

/ 

i-l 


exp^-v2/vs2j 


(C-9) 


Taking  the  density  moment  of  the  perturbed  distribution  function  and  sub¬ 
stituting  into  Poisson's  Equation, 
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k^1  *  -4ne(nle  -  n^)  ,  (C-10) 

results  in  the  following  dispersion  relation  for  electrostatic  inodes. 


0=1+ 


s  K  aDs  I  p=» 


!  ,  M  +  ivs  ~  ^Ds  z 

(2)*  K  vs 


/“  +  1vs  -  pn5\ 
\  (2)i  kz  vs  / 


W  exP(-xs) 


p=o 


1  + 


(2)*  K  vs 


ivs  >  +  ivs  -  pns\ 

'  (2)*  *,  ) 


!p(Xs)  exp(-A$)J 


-1 


where 


k  2  v  2 
>  _  Kv  vs 

'Ds  '  ’  s  "  2 

4-rrn  e  n 


2.  Ts 


Z(x)  =  -2  exp(-x2)  f  dt  expCt2)  +  exp(-x2)  . 

Jn 


(C-ll) 


(C-12) 


In  Eq.  (C-12),  Ip  is  the  modified  Bessel  function  of  the  first  kind  and 
order  p. 

The  incorporation  of  electron  collisions  with  ions  in  determination 
of  nle  and  of  ion  collisions  with  electrons  in  determination  of  n^.  poses  a 
special  problem  in  that  each  species  has  a  velocity  which  oscillates  as  part 
of  the  wave.  Hence,  the  collisional  contributions  for  each  species  on  the 
righ  hand  of  Eq.  (C-5)  should  depend  on  the  wave  fluid  velocity  of  the  other 
species.  However,  the  collisional  model  of  Eq.  (C-5)  does  not  include  this 
behavior. 

The  effect  of  interspecies  coulomb  collisions  on  drift-dissipative 
modes  can  be  estimated  from  the  ions  and  electron  fluid  equations  in  the 
limit  k^.  «  1.  For  simplicity,  we  take  the  zero-order  (background)  ion 
and  electron  velocities  to  be  zero.  Define  the  first-order  velocities  [vary¬ 
ing  as  expOk^y  +  ikzz  -  iwt)]  as  v+  and  v_  for  ions  and  electrons,  respec¬ 
tively.  The  first-order  fluid  velocity  equations  are 


9v+  q  E  VP. 

Tf  *  v±  x  %  ‘  V+  -  Vc±(v±  -  V  +  V  ‘  M 


(C-13) 


where  vc  denotes  the  collision  frequency  for  Ions  on  electrons,  and  v£ 
denotes  the  frequency  for  electrons  on  Ions.  Now  expand 

u  *  yO  +.  yl 
V+  Y+  + 


where 


0  -*■  3+  ^+- 

Vx  *  m — ~  + 

+  c  ±  m+  m+n 


v+  x  “c*  8  vn*l 


(\  X  “c.)  ♦  V  j(V?  -  V;  )  *  “«  +  *  "t 


(Effectively,  a  series  In  Inverse  powers  of  B  Is  generated  for  v+  ,  the 
component  of  wave  fluid  velocity  perpendicular  to  B.)  We  note 

V8l“cJB/|B|«c  8  -I«CJB/|B|  . 

After  some  algebra,  one  obtains  for  electrons 


q. 

vlx*  m 


V *  k  Vc-  +  !M 
v  -«  m-nl 

t„  Te  ikz  n,  t  q  II.,  ♦, 


m,.  nft  va  m_  v 

e  0  e  e  e 


and  for  Ions 


,  .  a* ,  x  k .  7p+x  *  v"+  f ■  Ei ,  7p+l  / vp*  pp.  \ 

*  *  c  c  L  ♦  +  J  >*  %.  ”+"“0+7 

»  (  "♦  E1  ,  pp+  \ 

\  *tt  ■+  "“c+  / 

Here,  k  s  B/|B|.  On  neglecting  Ion  velocity  ||  B  and  using  the  first-order 
electron  and  Ion  continuity  equations. 


♦  v  .  (n+  v+)  »  0 
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one  obtains  for  electrons 


*le  q*l 


ip  Icl  T 


(«0e  *  l6) 


b'e  (w  +  16  +  ID.  k  2 

*  V 


(C-14) 


and  for  Ions 


Uli  a  J!!L_  (u>De  “  p12  ~  un+  kf  pi? 

+  ky2  pi2  ♦  ID,  k  2) 


2  2  Te, 


n0  kb  Te 


(C-15) 


Here, 


°*  me  “e  3*  "o 


)a.  VT<  +  >>  »  V  kb<T1  ♦  V  , 

IV  T^J— m|»c  |-  ■*  *,e  “f  o  * 


VV 


qeB 

ft  = 

e  me 


6  = 


kbW 

m  v 


We  further  note 


el 


v-  ■  ven  +  ».i 

(C-15a) 

ven  *  3-5  *  ,0'9  % 

(C-15b) 

\  ■  -In  +  -3  -,1  k/  »,2 

(C-l 5c) 

v1n  ‘  8'3  *  M*'° 

(C-15d) 

ne[23.4  -1.15  log1Q  ne  +  3.45  log1n  TJ 

(C-15e) 

3  .  5  x  lO5'  T  3/Z 
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(C-15f) 


4 


T  3/2  ma 

v.  ~  nj  .  (=&)  — 

le  ei  'TV  m^ 

ne(mp/mj)^[23,4  -  1,15  lo9i0  ne  *  3.45  log1Q  Tg] 
2.12  x  107  T3/2  "" 


(C-15g) 


The  ratio  of  proton  to  ion  mass  is  written  as  nip/m.. .  The  quantities  vgn, 
vin*  vie'  vei*  and  vii  denote  the  electron-neutral ,  ion-neutral,  ion- 
electron,  electron-ion,  and  ion-ion  collision  frequencies,  respectively. 

In  formulas  (C-15a)  through  (C-15c),  all  quantities  are  in  cgs  units,  ex¬ 
cept  for  the  temperature,  which  Is  in  units  of  electron  volts.  The  form¬ 
ulas  (C-15a)  through  (C-15c)  are  taken  from  Ref.  10,  while  formulas  (C-15d) 
through  (C-15g)  are  taken  from  Ref.  31. 

Kinetically  from  Eqs.  (C-10)  and  (C-ll)  on  taking  p  =  0  contributions. 


In  the  limit  (w  +  ivi)/(2i  kz  v^)  »  1  on  neglecting  non-dominant  terms  in 
A..,  this  becomes 


If  one  identifies  ,  expands  to  first  order  in  a^,  and  notes 

wDi^RT  *  Ec,s*  and  (c-l6)+c,l^er  only  in  the  term  in  Dx  ky2  in  the 

denominator  of  (C-15).  Therefore,  we  modify  the  kinetic  result  (C-16)  with 
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7 

<*>-*■«+  ky  In  the  denominator,  so  that 


nii .  *1 


(id  +  1Vj)[l  ~  ^p(^^)e~  ]  "  ^  Ifl(*f)e 

<*>  +  iv^l  -  Iq(x1  )e  +  1D±  ky2 


(C-1 7 ) 


We  note 


=  -li1  I  VV!  j  >  *3  vii(ky  p1)2(ky  p1)2  .15(3.5  x  105)  (%*(%,  k  .5 

2  ~  vie(ky  p,)Z{1  +  T./Te)  (2.12  x  107)  'V  'V  y  Pl 


Dj-  ky 


Using  mp  =  16, 


r  >  18(ky  Pi)2  .  (C-18) 


Hence,  the  additional  term  is  non-dominant  for  ky  p1  >  .24,  and  it 
appears  reasonable  to  use  Eq.  (C-1 7)  for  all  ky  p^  <1. 

Taking  quasi -neutral ity  (nie»nH)  and  using  Eqs.  (C-14)  and  (C-17), 
with  a  *  IQ(x1 )  exp(-x..)  results  in  the  equation 


,  2 


<?-)  +  (br  +  1b,  H*-)  +  cr  +  1c,  -  0 


(C-1 9) 


where 


u'  =  id  +  10  k. 


b  _  +!De  n 
br“  +  v1  rt 


(C-20) 


cr  =  - 


n  k  i 

—  (1  +  rt,  — JL.  _L 
v1  T  Vl  v1 


1  Vi  '  \  1-0  V]  J 
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This  equation  can  be  solved  by  the  standard  quadratic  formula  to  yield  two 
values  for  which  by  Eq.  (C-20)  determine  values  for  u>.  Our  computer 
program  calculates  the  root  with  largest  algebraic  Imaginary  part  (=y).  If 
y  >  0,  a  diffusion  coefficient  Dx  *  y/kx2  with  =  (p^  Lx);  If  y  <  0, 
no  diffusion  Is  anticipated. 

For  a  given  set  of  background  physical  parameters  (e.g.,  ambient 
magnetic  field.  Ion  species.  Ion  temperature,  electron  temperature,  elect¬ 
ron  density,  neutral  density,  background  plasma  scale  length  parallel  to  the 
ambient  magnetic  field,  plasma  gradient  scale  length  perpendicular  to  the 
ambient  magnetic  field),  falsest*  ky  are  scanned  In  geometric  Increments  over 
the  range 

“cl  1 

rr1  >  k  >  X  , 


and  values  of  kz  are  scanned  in  geometric  Increments  over  the  range 


min  (-S-  .  — )  >  K  >  r1 


Further,  It  Is  necessary  to  check  >  w. 

The  larges  value  of  y  for  a  pair  of  ky,kz  over  the  permissible  range 
°f  ky,kz  parameter  space  Is  taken  to  determine  the  diffusion  coefficient. 

Some  further  insight  into  the  effect  of  Interspecies  coulomb  inter¬ 
actions  is  provided  by  consideration  of  the  limit  ky  »  1  »  ky  pg.  Then 
one  expects  from  finite  ion  Larmor  radius  effects  that  |v+|«|v_|  and 

v  *  v  .  +  v  in  Eq.  (C-13).  Also,  the  Ion  behavior  should  be  predomin- 
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antly  Independent  of  v^e  because  of  the  dominance  of  ion-ion  viscosity  over 
Ion-electron  collisions  (see  Eq.  (0-18)). 

On  using  Eq.  (C-13)  for  the  electrons,  we  have 


+  16  +  i(\ 


+  vc-)(ky  pe^ 


u  +  1«  +  1{vn_  +  vc_)(ky  Pe)‘ 


(C-21) 


while  for  the  lor',  we  have  Eq.  (C-16). 

On  taking  quasi -neutrality  and  using  Eqs.  (C-16)  and  (C-21),  the  solu¬ 
tion  equation  Is  again  of  the  form  of  Eq.  (C-19)  with  6  replaced  by 

6  +  (v  +  v,  )(kupj2  and  0,  replaced  by  0  in  Eq.  (C-20). 
n-  c-  y  e  j. 


In  this  case,  electron-ion  collisions'!  effects  are  neither  Inherently 

stabilizing  nor  destabilizing.  For  kz  such  that  k  2  kb  Tg/ine  vg  «  |w *|, 

the  effects  are  destabilizing;  for  k  such  that  k  *  k.  T  /m  v  %>  I, 

z  7  «  «  e  e  1  De  1 


the  effects  are  stabilizing. 


APPENDIX  D 


SURVEY  OF  OTHER  DENSITY-GRADIENT  INSTABILITIES 
POSSIBLY  RELEVANT  TO  ANOMALOUS  DIFFUSION  UNDER  HANE  CONDITIONS 


Anomalous  diffusion  In  plasmas  is  not  a  priori  limited  to  diffusive 
behavior  driven  by  any  particular  plasma  mode.  Hence  It  is  relevant  to  con¬ 
sider  modes  In  addition  to  the  drift-dissipative  mode^'^s  possible  sources 
of  wave  turbulence  and  diffusion. 

Our  purpose  here  Is  to  survey  the  various  types  of  linear  wave  insta¬ 
bilities  (besides  the  drift-dissipative  mode)  which  could  be  driven  unstable 
as  a  result  of  the  free  energy  associated  with  an  ambient  density  gradient 
and  with  both  an  ambient  density  gradient  and  an  ambient  electric  field  for 
lower-hybrid  waves.  We  consider  Ionospheric  situations  with  the  following 
parameters  relevant  to  long-time  strlation  behavior  and  scintillations: 

6  -3 

n  *  n.  =  10  cm  ,  T  *  T*  ■  0.1  ev, 
e  i  e  i 

B  *  0.3  gauss  •  (D-l) 

In  Equation  (D-l)  the  quantities  n,  T,  and  B  denote  the  density, 
temperature,  and  ambient  magnetic  fields,  respectively.  The  subscripts  "e" 
and  "1"  represent  electrons  and  Ions,  respectively,  a  quasineutral  plasma 
made  up  of  singly  Ionized  oxygen  Is  also  assumed.  Given  the  parameters  of 
Equation  (D-l)  the  ratio  of  plasma  to  magnetic  field  pressure,  B,  Is 

B  =  9  x  10”5  (D-2) 

which  Is  comparable  to  the  ratio  of  electron  to  oxygen  masses, 

■jp  ■  3.4  x  10'5  .  (D-3 ) 

m1 

Consequently,  except  for  some  possible  modes  (to  be  discussed  subsequently) 
with  wavelengths  parallel  to  the  magnetic  field  of  the  order  of  several  tens 
or  hundreds  of  kilometers,  the  relevant  modes  from  the  standpoint  of  anomalous 
transport  will  primarily  be  electrostatic. 
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We  assume  a  z-directed  magnetic  field  of  strength,  B,  an  x-directed 
density  gradient  and  uniform  plasma  properties  in  the  y-di recti  on.  This 
slab  like  plasma  model  most  closely  approximates  the  true  ionospheric  plasma 
if  the  wavelengths  of  modes  are  much  smaller  than  the  inhomogeneity  scale 
lengths.  This  geometry  is  shown  in  Figure  30. 

In  the  following  subsections  we  analyze  the  linear  stability  of  drift 
instabilities  associated  with  cyclotron,  acoustic,  lower-hybrid  and  electro¬ 
magnetic  behavior.  The  calculations  support  the  conclusion  that  the  drift- 
dissipative  mode  is  the  dominant  wave  turbulence  factor  in  anomalous  diffusion 
under  HANE  conditions. 


A.  ANALYSIS  OF  MODES 


1.  Ion-Cyclotron  Drift  Modes 

Drift  instabilities  with  frequencies  comparable  to  the  fundamental  or 
harmonics  of  the  ion-cyclotron  frequency  can  be  driven  unstable  if 

P^i  <  |wD|  (D-4) 

where  "p"  is  equal  to  a  nonnegative  integer,  a>D  is  the  diamagnetic  drift  fre¬ 
quency,  e.g., 

“D  '  VlZ/C1L„ 


with  v^  *  T^/m.j  ,  n.  =  eB/m^c  . 

(5) 

An  equivalent  condition  to  Equation  (D-4)  is 

p  <  (fcwpJ  r~  (D-5) 

y  n 

Hence,  for  situations  with  weak  density  gradients  such  that 

Pi  «  Ln  (D-6) 

we  require 

kyPj  »  1  (D-7) 

for  instability. 
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A  factor  which  can  severely  limit  the  destabilization  of  ion-cyclotron 
modes  is  ion-ion  and  ion-neutral  collisions.  In  particular  if (32, 33) 

Q1<v1  5  vin  +  -3  k22  pi2  vii  (D-8) 

ion-cyclotron  modes  are  no  longer  effective  in  the  plasma  dynamics.  In 
Equation  (0-8)  we  have  inserted  the  .3  in  the  ion-ion  contribution  to  the 
effective  collision  frequency  as  a  reflection  of  the  results  of  a  detailed 
calculation  using  the  Fokker-Planck  operator  in  the  k2p.«l  limit/11  1 

Using  Equations  (D-5)  and  (D-8)  for  singly  ionized  oxygen  ions  with 

nn  =  108  cm"3  ,  ng  =  106  cm'3  ,  Tg  =  .1  ev  , 

B  =  .3  G  (0-9) 


it  is  apparent  that 

kj_Pi  <  10  ,  Pi/Ln>  .1  (D-10 ) 

is  necessary  for  the  existence  of  ion-cyclotron  modes.  Hence,  the  ion-cyclotron 
drift  instabilities  require  rather  sharp  density  gradients. 

The  conditions  given  by  Equation  (D-10)  represent  necessary  criteria  for 
the  existence  of  the  ion-cyclotron  instabilities.  However,  even  when  inequality 
(D-8)  is  not  satisfied,  ion  collisions  can  still  stabilize  ion-cyclotron  modes. 
Using  the  general  dispersion  relation  for  electrostatic  modes  given  by  Eq.  (C-ll) 
in  Appendix  C  and  assuming  short  perpendicular  wavelength  perturbations  such 
that 

kJ_pi  ^  1  ( D - 1 1 ) 

we  can  rewrite  the  dispersion  relation  as 


0 


1 

K  ui 


w  +  iv..  -  pfii 

i - 

<*>+  iv^ 


(2t r)h  k±p.  J 

+  lv.  1 

-  p (2tt)^  k^p. 


(D-12) 
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where  pfl^  Is  the  harmonic  of  the  Ion-cyclotron  frequency  to  which  w  is  most 
closely  approximated  (e.g.,  and  xe  Is  the  electron  suscepti¬ 
bility.  From  Equation  (D-12)  the  expression  for  the  frequency  Is  given  by 


where 


w  *  pfi.  + 


*01  +  *1 


1  {2^  kj_Pl 


xer  +  2"  "  *xei 

k  ADi 


-  r 


(*er  +  rr  ?)  +  xei 

k  ad1 


7T1" 

K  Aq1 


(D-13) 


iuj. 


Xer  3  ffe  (xe)  .  Xet  3  Im 


Ion  collisions  tend  to  stabilize  the  ion  cyclotron.  For  situations  when 
Instability  can  occur, 

lwDlf  >  IP°f  I  •  (0-14) 


Stability  is  assured  If 


2  2 

.  3kj^  p.|  v  ^  j 


2u)i 


01 


(2nP  k 


clp1 


=  (?) 


P1 


(D-15) 


Using  Equation  (D-3)  and  assuming  the  parameters  given  by  Equation  (D-9)  as 
well  as 

P1/Ln  3  -1  (D-16) 

we  see  that 

kj.Pt  >  3  ( D- 17) 

for  stabilization  by  Ions.  Even  for  electron  densities  as  low  as  10®  cm-3 
the  requirement  that 

k  Pt  >  9  (D-18) 


for  stability  shows  that,  unless  density  gradients  are  very  sharp  or  densities 
are  very  low.  Ion-cyclotron  drift  modes  are  not  likely  to  be  driven  unstable 
In  the  Ionosphere  solely  by  density  gradients. 
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2.  Acoustic  Modes  (u)>fl.) 

Acoustic  waves  with  frequencies  above  the  Ion-cyclotron  frequency 
satisfy  the  following  expression  for 

\*e  >>  kzve:  (D-19) 


The  criterion  for  the  derivation  of  Equation  (D-20)  from  Equation  (C-ll)  is  that 

kjPj  <  p  (D-21) 

where 

a)  s  pQ.  .  (0-22) 

In  the  Ionosphere  the  electron  and  ion  temperatures  are  comparable  and  we 
requi re 

Pi  >  Ln  ( D-23 ) 

for  instability.  It  then  follows  that  the  ion  Larmor  radius  must  be  larger 
than  the  density  gradient  scale  length, 

>  a)  ( D— 24 ) 

for  a  gradient  driven  acoustic  instability.  Hence,  the  acoustic  instability 
will  occur  only  for  situations  where  the  density  gradient  is  extreme.  Since 
only  ion  dynamics  have  been  used  to  establish  Equation  (D-24),  our  arguments 
also  apply  to  the  case  when  Equation  (D-19)  is  not  satisfied  and  the  electrons 
are  collisionless. 

3.  Lower- Hybrid  Waves 

In  the  coll  is  ionless  limit  the  dispersion  relation  for  lower-hybrid 
waves  satisfying  Equations  (D-21)  and  (D-22)  is 
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Again,  as  for  the  case  of  acoustic  modes,  we  require  that 

P1  >  Ln  (D-26) 

be  satisfied  for  instability.  Hence,  as  in  the  case  of  acoustic  modes, 
collisionless  lower-hybrid  waves  driven  solely  by  the  density  gradient  and 
satisfying  Equations  (D-21)and  (0-22)  are  not  likely  to  be  present  in  the 
ionosphere. 

The  case  of  a  background  electric  field  has  recently  been  raised  by 
(341 

Huba  et  al.v  '  in  connection  with  the  interpretation  of  enhanced  backscatter 
from  Equatorial  Spread  F  by  the  ALTAIR  radar  on  Kwajalein.  The  calculation 
was  done  assuming  that  the  inequality  represented  by  Equation  (D-21)  was 
reversed  (i.e.,  k^p^  >  p  for  w  a  pfi. ).  Under  such  circumstances  the  quantity 
p^/Ln  can  be  less  than  one  for  a  lower-hybrid  instability.  However,  the 
stabilizing  influence  of  electron  collisions  discussed  here  are  not  included 
in  Reference  (1)  and  so  the  examples  used  in  the  reference  to  explain  the 
backscatter  are  Inappropriate. 

The  plasma  model  of  Reference  34  assumes  a  uniform  z-dlrected  magnetic 
field  of  strength  B,  uniform  plasma  properties  in  the  y-direction,  and  an 
x-directed  density  profile  of  the  form 

n0(x)  =  n0(l-x/Ln)  .  (D-27) 


The  ions  are  assumed  to  be  In  local  thermodynamic  equilibrium  with  an  ambient 
electric  field  of  strength,  Eq(x)  and  potential  cpQ (x) : 

iLvm  (n0(x)>  =  E0(x)  *  -  V*0(x)  (D-28) 

where  is  the  ion  temperature  in  electron  volts. 


Ill 


To  lowest  order,  the  ion  dynamics  can  be  treated  as  being  straight- 
line  and  unmagnetized  if  the  ion-ion  collisional  viscosity  decorrelates  the 
ion  gyromotion  over  a  scale  length  on  the  order  of  an  inverse  y-directed 
wave  number,  k  Quantitatively  for  the  ions  to  be  unmagnetized, 

•3(vii/fli)  ky2pi2  >>  1  (D-29) 

where 

vii  =  f  ^  ^3/1  *-23’4  '  1,15  1og10^no^  +  3,45  1og10^Te^  ’ 


n  --SL 

“i  m^c  ’ 


p<*7 


vi 


(D-30) 


are  the  ion- ion  collision  frequency,  ion  gyrofrequency,  ion  gyroradius,  and 
ion  thermal  speed,  respectively.  In  the  definitions  of  Equation  (D-30),  the 
ions  have  been  assumed  to  be  singly  ionized. 


Contrary  to  the  ion  motion,  the  gyration  of  the  electrons  in  the  ambient 
magnetic  field  is  crucial  to  the  development  of  a  lower-hybrid  instability  and 
is  necessarily  included  in  all  treatments  of  lower-hybrid  waves.  What  are  not 
included  in  the  treatment  of  Reference  34  are  the  electron-electron  and  elec¬ 
tron-ion  collisions  which  obviously  must  be  considered  in  a  self-consistent 
collisional  treatment  of  lower-hybrid  waves  in  the  ionosphere  if  ion-ion  col¬ 
lisions  are  included.  The  electron-electron  and  electron-ion  collision  fre¬ 
quencies  (va_  and  v  respectively)  have  the  form 

CC  v  • 


v 


ee  ”  vei 


vii 


(D-31) 


For  situations  where  Equation  (D-29)  is  satisfied  and  the  ion-ion  mean  free 
path  is  much  larger  than  the  inverse  y-wavenumber. 
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(0-32) 


r 


v1/v1i  ~  ky 


the  ion  dynamics  in  the  lower-hybrid  wave  are  collisionless  and  are 
described  by  the  1 inearl ized  and  unmagnetized  Vlasov  equation 


ar 


e  h  a  .  /-*•-*•* 
"“T  '  3?  f0'  r,v 


(D-33) 


In  Equation  (0-33),  wave  and  background  quantities  are  represented  by 
subscript  ”1"  and  ”0",  respectively,  and 


fjj(r»v»t)  *  f^(v)  expCk^y  -  ut)] 


fQi(r,v)  =  (2ir  T./m^.)"372 


(D-34) 


V 


-y  i  k  *  exp[i(kyy  -  mt)] 


The  electron  dynamics  in  the  wave  are  magnetized  and  collisional. 
For  lower-hybrid  waves,  the  electron  collisions  with  ions  can  be  assumed 
to  be  Lorentz-like,  and  the  effect  of  electron-electron  viscosity  can  be 
assumed  to  be  additive  to  the  electron-ion  collision  frequency  in  the 
following  way: 


where 


ve  =  vei  +  ,3vee  {ky  pe}‘ 


(0-35) 


(D-36) 


In  Equation  (D-35),  the  effective  electron  collision  frequency  is 
defined  by  vfi.  An  appropriate  linearized  kinetic  equation  for  the  elec¬ 
trons  is 
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Here, 


nle(r,t)  -  n]e  exp[1  (ly  -  u»t)3 
f1e(?,v,t)  *  fletf)  expfltkyy  -  wt)] 


^(r.v)  *  ( 2-jt  TA)3/2  exp|- 


e'  e' 


me  v‘ 


In  the  Ionosphere,  the  relation 


2  2  2 
<*»pe  4irn0  e  /Ne  4™Q  c  me 

^  2  ~  rt  /_  _  \  2  _2 


»  1 


ne  (e  B/mecV 


B 


is  satisfied.  Since  for  the  lower-hybrid  mode  one  also  has 

|ne|  »  to 

the  mode  can  be  assumed  to  be  quasi -neutral . 

Using  standard  techniques,  the  dispersion  relation  for  the  lower- 
hybrid  modes  can  be  written,  using  Equations  (D-32)  and  (D-37),  as 


0  =  [1  +  yZ(y)] + 


f,  H 

u)  +  1 


L  *  ~  “De  + 
[  “  +  iv£ 

i-  ■ 
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Here, 


Z(y)  =  -2  exp(-y2) /  dt  exp(t2)  +  if*)*  exp(-y2) 
0 


Co  -  ta)i 


01 


<2>*  ky  V1 

Ti  ,  Te 

W01  *  ky  mi  n.  Ln  ’  wDe  =  ky  mfi  Ln 


(D-38) 


(D-40) 
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r0  =  Vxe*  exp^"xe^  »  xe  =  ky2pe2  *  0>-41) 

and  IQ  is  the  modified  Bessel  function  of  the  first  kind  and  zero  order.  In 
writing  Equation  (0-40),  Eq.  (D-28)  and  a  frame  of  reference  moving  with  the 
ExB  velocity  have  been  used. 

Assuming  all  the  ions  with  the  exception  of  resonant  ions  are  in  local 
thermodynamic  equilibrium  with  the  wave, 

!<*>  -  uiDi!  «  (2)1  kvi  ,  (D-42) 

the  dispersion  relation.  Equation  (D-40),  can  be  rewritten  as 


0  »  w(2  -  rQ)  -  a>Di  +  1(J)*  a,  +  2i  ve(l  -  x0)  •  (0-43) 

If  in  addition  the  effective  electron  collision  frequency  and  the  growth 
rate,  y,  are  much  smaller  that  the  real  part  of  the  frequency,  w, 

lY.ve(1-r#)|  «  M  ,  (0-44) 


then 
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(0-45) 


The  last  term  on  the  right-hand  side  of  Equation  (D-45)  represents 
the  stabilizing  influence  of  electron  collisions  on  the  lower-hybrid  drift 
wave  and  was  not  included  in  the  treatment  of  Reference  34.  Hence,  unlike 
the  case  of  Reference  34,  the  consistent  treatment  of  collisions  results 
In  a  threshold  criterion  for  the  density  gradient  scale  length  necessary 
for  instability.  This  criterion  is 


( 


2  me^  2  ~  r0 

*"1  (|kpe|r0)* 


(D-46) 


Example 

As  an  example,  consider  an  ionospheric  plasma  with 


% 

B  =  .3G  ,  ir  =  3.4  x  10 

m1 


106  cm'3 


vee  =  1180  Hz  *  Te  =  Ti  =  -1  eV  • 
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T 


The  criterion  of  Equation  (0-46)  Implies  that: 

p-  >  0.56  for  |k  p J  *  .1  , 

Ln  y  e 

A  >  0.37  for  | k  p  I  =  .7  , 
n  J 

Pi 

L1  >  0.44  for  jk  pe|  =  1.0  .  (D-48) 

n  y 

The  values  of  the  y-wavenumber  such  that  |k  p  j  =  .1  and  .7  are  appropriate 
to  radar  backscatter  at  100  MHz  and  50  MHz,  respectively. 

From  the  criterion  of  Equation  (D-40),  it  is  apparent  that  the  examples 
of  Reference  34  with  the  same  parameters  as  those  used  here  but  with  p./Ln  = 
0.026  actually  correspond  to  a  stable  lower-hybrid  drift  wave.  Based  on  Equa¬ 
tion  (D-46),  it  is  clear  that  for  the  lower-hybrid  drift  instability  to  be 
driven  unstable,  the  plasma  density  would  have  to  be  reduced  to  ^10^  cm-3  if 
other  plasma  parameters  were  unchanged.  For  such  densities,  Equation  (D-29) 

would  still  be  satisfied  as  long  as 

|ky  Pel  »  0.1  .  (D-49) 

Hence,  it  appears  that  lower-hybrid  drift  modes  are  not  likely  to  be 

A  o 

sources  of  anomalous  diffusion  under  HANE  conditions  (nQ  »  1(T  cm" J),  unless 
gradient  scale  lengths  are  of  the  order  of  10  m  or  less. 

4.  Electron-Cyclotron  Modes 

Electron-cyclotron  waves  have  frequencies  comparable  to  the  funda¬ 
mental  or  a  harmonic  of  the  electron-cyclotron  frequency.  Analogous  to  the 
treatment  for  ion-cyclotron  drift  instabilities,  we  require  that 

kx  Pe2/ln  >  P  (D— 50) 

for  instability.  In  Equation  (D-50),  the  quantity  p  is  the  electron- 

cyclotron  harmonic  number  to  which  the  wave  frequency  is  most  closely 

approximated.  Again,  as  in  the  case  of  ion-cyclotron  modes,  we  require 

* 

•  3  kj2  peZ  vee  «  |fiel  (0-51) 
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for  the  existence  of  electro-cyclotron  modes.  In  terms  of  Ion  quantities, 
and  using  Equation  (D-50),  we  find  that  a  necessary  condition  for  the 
destabilization  of  electron-cyclotron  drift  waves  Is 

.3  kj  p*  Vi.  <  ni(me/2m.)i  .  (D-52 

Hence,  the  density  constraints  for  the  existence  of  electron-cyclotron 
drift  Instabilities  are  more  stringent  than  those  for  ion-cyclotron 
Instabilities  by  the  ratio  (me/2m^. 

5.  Electron-Acoustic  Modes 

The  destabilization  of  electron-acoustic  modes  with  frequencies  well 
above  the  electron-cyclotron  frequency  would  require 

|kxPe!«p  (D-53 

where  p  is  an  Integer  such  that 

to  =  Pfig  .  (D-54 

Hence,  following  the  arguments  used  for  ion-acoustic  modes,  it  follows 
that 

M  >  Ln  *  (D-55 

which  corresponds  to  a  very  extreme  condition  on  the  density  gradient. 

Hence,  electron-acoustic  modes  driven  unstable  by  density  gradients  should 
not  exist  in  the  atmosphere., 

6.  Electromagnetic  Drift  Modes 

As  the  plasma  energy  becomes  larger  in  a  plasma  configuration  with 
a  fixed  magnetic  field  strength.  It  becomes  easier  for  the  plasma  to  modify 


the  ambient  magnetic  field.  Under  such  circumstances,  instabilities  in  a 
plasma  can  have  electromagnetic  properties. 

We  have  shown  that  electrostatic  modes  driven  unstable  solely  by 
density  gradients  In  the  Ionosphere  are  generally  expected  to  have  frequencies 
which  are  less  than  the  Ion-cyclotron  frequency.  Similarly,  because  the  wave 
frequency  must  be  smaller  than  the  diamagnetic  drift  frequency  for  instability 


any  possible  electromagnetic  modes  driven  by  density  gradients  in  the  ion¬ 
osphere  should  also  have  frequencies  below  the  ion-cyclotron  frequency. 

Electromagnetic  modes  with  frequencies  below  the  ion-cyclotron 
frequency  and  driven  unstable  by  a  density  gradient  have  the  property  that^ 


*  VA  '  T'  Ln 


(0-56) 


where 


8*  n0s  Te 


(4tt  nQ  mf)‘ 


(D-57) 


are  the  ratio  of  the  species  s  pressure  to  magnetic  field  pressure  and  the 
Alfyen  speed,  respectively.  The  quantity  $  is  ubiquitous  in  discussions  of 
low-frequency  electromagnetic  instabilities.  The  low-frequency  electromag¬ 
netic  instabilities  also  have  the  property  that 


i  -  — 


k2  c2 


(D-58) 


where  e„jtl  is  derived  in  Appendix  E  and  the  speed  of  light  is  written  as 
c.  Using  Equation  (E.4),  we  can  verify  Equation  (D-58)  if 


Ui 

>-l*  > 


k  2  c2 


( D-59) 


Substituting  the  diamagnetic  drift  frequency  for  the  frequency  in  Equation 
(0-59)  and  using  Equation  (D-56)  allows  for  a  criterion  on  the  parallel 
current  to  be  written  for  the  electromagnetic  modes  to  be  written 


For  the  parameters 


B  «  .36 


X,.  >  2tt 
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pe  1 
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(D-61 ) 
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In  the  Ionosphere  with  singly  ionized  oxygen. 


A 


II 


(D-62) 


Unless  the  quantity,  Lp/P^.  is  much  smaller  than  10,  the  necessary  requirement 

on  the  parallel  wavelength  for  electromagnetic  instabilities  seems  prohibitive. 

8  -3 

Even  for  the  ion  and  electron  densities  of  10  cm  ,  we  require  the  seemingly 
difficult  inequality 


to  be  satisfied. 


(0-63) 


However,  it  is  interesting  to  note  that  for  plasma  parameters  such  as 


B  *  .03G 


"0e  n0i  ‘  '°8  Cl"’3 


.1  eV  , 


ee  *  Pi  . 

there  is  the  possibility  that  electromagnetic  modes  will  be  driven  unstable 
by  the  density  gradient.  Necessary  conditions  for  instability  are 

Be  *  8i  85  .45  ,  At|  >  km  •  (0-65) 

In  this  case,  the  total  plasma  energy  is  only  slightly  smaller  than  the  mag¬ 
netic  field  energy: 

e  =  Be  +  Bi  *  .90  .  (D-66) 


B.  SUMMARY 

It  is  possible  that  lower-hybrid  and/or  electromagnetic  drift  modes  can 
be  driven  unstable  under  ionospheric  conditions.  However,  in  both  cases  there 
appear  to  be  rather  stringent  limitations  on  the  parametric  range  of  instability. 
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APPENDIX  E 

DERIVATION  OF  THE  CURRENT  PRODUCED  IN  A  WAVE  WITH  u  «  Q. 
IN  THE  DIRECTION  OF  THE  AMBIENT  MAGNETIC  FIELD 


In  this  appendix,  we  derive  the  current  produced  by  a  low-frequency 
wave  (i.e.,  w  «  in  the  direction  of  the  ambient  magnetic  field  (which  is 
taken  to  be  the  z-direction) .  We  assume  a  collisional  plasma  (i.e.,  ve.  »  U) 
and  use  the  fluid  equations. 


Here, 


0  =  -ik2  Tq  n,  -  e  nQ  Ezl  -  v0  nQ  vzl  , 

n  cE  ,  3nQ 

0  =  "1  nl  +  "B^  17  +  lkz  n0  vzl 


V  =  V  +  v  • 

e  en  ei 


solving  for  the  perturbed  velocity  parallel  to  the  magnetic  field,  we 
obtain 


fzl  ‘  '  m 


■  ^  [.  * «,  '.fe  ■  ”,  i^r-) 

f  - 


The  parallel -parallel  component  of  the  dielectric  tensor  can  then  be 
written  as 


(E-la) 


(E-lb) 


(E-2) 


(E-3) 
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In  Eq.  (E-4), 
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